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Abstract 

Simple optical instruments are linear optical networks where the incident 
light modes are turned into equal numbers of outgoing modes by linear trans- 
formations. For example, such instruments are beam splitters, multiports, 
interferometers, fibre couplers, polarizers, gravitational lenses, parametric am- 
plifiers, phase-conjugating mirrors and also black holes. The article develops 
the quantum theory of simple optical instruments and applies the theory to 
a few characteristic situations, to the splitting and interference of photons 
and to the manifestation of Einstein-Podolsky-Rosen correlations in paramet- 
ric downconversion. How to model irreversible devices such as absorbers and 
amplifiers is also shown. Finally, the article develops the theory of Hawking 
radiation for a simple optical black hole. The paper is intended as a primer, 
as a nearly self-consistent tutorial. The reader should be familiar with basic 
quantum mechanics and statistics, and perhaps with optics and some ele- 
mentary field theory. The quantum theory of light in dielectrics serves as 
the starting point and, in the concluding section, as a guide to understand 
quantum black holes. 



1 Introduction 



Consider a semi-transparent mirror, the glass of your window, for example. The 
mirror partially reflects light and is partially transparent. If the material of the 
mirror is not absorptive the incident light is exactly split into the reflected and the 
transmitted component. Now, light consists of photons, of indivisible light particles. 
How does this beam splitter act on individual photons [32* | I14fi| ? How are photons 
split? Or, in another experiment [71], suppose you take a semi-silvered mirror with 
50:50 transmission-reflection ratio. You let exactly one photon propagate towards 
the frontside of the mirror and you send another single photon towards the backside, 
and let them interfere. The interference between two light beams depends on their 
relative phase. If the phase difference is right, the sum of the two incident beams, 
the two photons, emerge behind the mirror and if they interfere with the opposite 
phase they appear in front of the mirror. But single photons are not supposed to 
carry a precise phase, because phase is a wave property and individual photons are 
particles. So what happens [TT j I146j ? 

Such conundrums are beginning to occupy people's minds for other reasons than 
purely academic curiosity, because they may fundamentally alter our approach to 
secure data communication EHj. Suppose that Alice wants to send a secret 
message to Bob, carried by photons in a glass fiber |174j or through space jHOj- Eve, 
the eavesdropper, tries to intercept the message without getting caught. Clearly, in 
order to do so, she must probe the stream of messenger photons, for example using 
something like a beam splitter. However, knowing the principal quantum effects 
of beam splitters, Alice and Bob may infer from the statistical properties of the 
transmitted photons that their secrecy is at risk and may discard the communication 
channel. In a more sophisticated eavesdropping attempt, Eve might amplify the 
incident light and extract, by beam splitting, the bits that are sufficient for her. 
The rest, with the same amplitude as the original, is transmitted to Bob. Would 
Alice and Bob notice Eve's subtle interception? 

The quantum physics of simple optical instruments, such as beam splitters and 
amplifiers, is clearly important when the quantum nature of light is used for practi- 
cal (or academic) purposes. Equally importantly, some seemingly simple questions 
about light and the relatively simple experiments to demonstrate their intriguing 
answers do both illuminate and challenge our understanding of the quantum world 
146J. Furthermore, in studying the quantum physics of simple optical instruments 
we may see connections to a much wider and sometimes quite exotic range of physics. 
For example, the Hawking radiation of black holes jHZI is related to the quantum 
optics of moving media and combines aspects of beam splitters and amplifiers. 

In this article we analyze the principles of quantum-optical networks where two or 
more beams of light interact with each other. The networks are assumed to be linear 
in the sense that the output amplitudes depend on the input amplitudes by a linear 
transformation, although the physics of such networks may be based on non-linear 
optics |2o1 1148| |168j . The article is intended as a primer, as a nearly self-consistent 
tutorial, rather than a literature survey. The reader should be familiar with basic 
quantum mechanics and statistics, and perhaps with optics and some elementary 
field theory. However, when appropriate, we quote the major mathematical results 
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needed, instead of deriving them, for not overburdening this article. We use the 
notation of the book (TlOj and some of the basic results of quantum optics explained 
there. First we begin with an example, quantum light in planar dielectrics. Then we 
extend the central features of this example to quantum-optical networks in general. 
We describe the quantum physics of such networks in the Heisenberg picture and 
in the Schrodinger picture, and with the help of quasiprobability distributions such 
as the Wigner function In Sec. 4 we develop the quantum optics of the beam 

splitter, because this simple device is the archetype of all passive optical instruments, 
and because the beam splitter is capable of demonstrating many interesting aspects 
of the wave-particle dualism. In Sec. 5 we analyze absorbers and amplifiers, which 
are irreversible devices, and show how they can be described by effective models, 
such as fictitious beam splitters and parametric amplifiers. In Sec. 6 we develop the 
quantum theory of parametric amplifiers and phase- conjugating mirrors. Parametric 
amplifiers have been widely used to experimentally demonstrate the nonlocality of 
quantum mechanics in versions of the Einstein-Podolsky- Rosen paradox [2111 EH]- 
Section 7 returns to the starting point, to quantum light in dielectric media. We 
show how moving media may establish analogs of black holes and we analyze the 
essentials of Hawking radiation |HZj . Throughout this article, whenever possible, we 
try to use models that are simple but not too simple. 

2 Quantum optics in dielectrics 

Many passive optical instruments such as lenses or beam splitters consist of dielec- 
tric materials like glass that influence the propagation of light without causing much 
absorption. Dielectrics are linear-response media — their effect on light is propor- 
tional to the electromagnetic field strengths. The field induces microscopic dipoles 
in the atoms constituting the dielectric medium. The dipoles constitute macro- 
scopic electric polarizations and magnetizations that are proportional to the applied 
electromagnetic field and which act back onto the field. An isotropic medium is 
characterized by two spatially dependent proportionality factors, the electric per- 
mittivity e and the magnetic permeability \x j74~ | 1100) . We assume that the medium 
is not dispersive and not dissipative. In this case both e and \i are real and do not 
depend on the frequency of light within the frequency window we are considering. 
The square root of the product of e and \i gives the refractive index that describes 
the degree to which the phase velocity of light in the medium deviates from the 
speed of light in vacuum, c. 

The quantum theory of light in dielectrics has been subject to a substantial 
literature summarized to some extent in Refs. [TS1 EH E3 E3 1185] . Traditional 
quantum optics is the subject of the recent books HS| H7J EDI 1123 1123 1131 
I149| I161| I167| I185| 1189] . Here we consider the simplest possible case where the 
dielectric functions e and \x vary in one direction of space only, say in x direction. 
Furthermore, we assume that the electromagnetic waves propagate in this direction 
as well and we select one of the two polarizations of light. In this way we arrive at 
an effectively one- dimensional model. 
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2.1 Classical fields 



Consider the classical electromagnetic field characterized by the electric field strength 
E and by the magnetic B field in SI units. The electromagnetic field obeys the Prin- 
ciple of Least Action [HZ] with the Lagrangian density [SI] 

In one spatial dimension the vector potential in Coulomb gauge [17] is effectively a 
scalar field A with 

E = -d t A, B = d x A. (2.2) 

Throughout this article we abbreviate partial-differentiation operators such as d/dt 
and d/dx by d t and d x . We obtain from the Lagrangian (J2.1j) the Euler-Lagrange 
equation 

-d x -d x A- — 

e fj, cr 

the wave equation of light in one- dimensional media at rest. We define the scalar 
product between two fields with vector potentials A\ and A 2 as 



d x -d x A--%A = 0, (2.3) 



(A 1 ,A 2 ) = 1 -^ J (Ald t A 2 -A 2 d t Al)edx. (2.4) 

As usual, h denotes Planck's constant divided by 2tt. The scalar product ()2.4|) is a 
conserved quantity as a consequence of the wave equation (|2.3j) . 

d t (Ai, A 2 ) = , (2.5) 

and the product plays an important role in the mode decomposition of quantum 
light. 



2.2 Quantum fields 

According to the quantum theory of light [4*7] I123| I127j the vector potential is 
regarded as a quantum observable, as a Hermitian operator A that depends on 
space and time. We represent A as a superposition of modes 

A(x, t) = J2 (Mx, t) a k + A* k (x, t) 4) . (2.6) 

k 

The mode functions A k satisfy the classical wave equation (j2.3j) . They describe 
how single light quanta, photons, propagate in space and time, given the initial 
and boundary conditions that determine the particular A%. They also describe how 
coherent states gll(J| I123| I127j propagate, states that describe classical light fields. 
The mode functions characterize the classical, wave-like, properties of light, whereas 
the quantum amplitudes ak describe the quantum features of light. We require that 
the mode functions are orthonormal with respect to the scalar product ([2.4)1 . 

(A fc ,Afc/)=<W, (Al,A k ,) = 0. (2.7) 
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As a consequence, the commutator relation between A and the canonically conjugate 
momentum |191j implies 

[ofe, a\] = 5 kk , , [a k , a k ,] = . (2.8) 

Therefore, light quanta are bosons [191 , i.e. quanta of harmonic oscillators, with 
annihilation operators a k and creation operators at. Each mode of light represents 
an electromagnetic oscillator, the best harmonic oscillators currently known (with 
the smallest anharmonicity, generated by the vacuum polarization due to electron- 
positron pairs [ZOl 1130] . an effect beyond our model). Furthermore, if we choose 
monochromatic mode functions with frequencies u> k , we obtain for the field energy 



/ 1 + f^ 2 ) dx = £ hUk ( a i a * + 1) 



(2.9) 



Each mode contributes to the total energy as the energy quantum ftw k times the 
photon number a\a k . The additional vacuum energy J2 k does not depend on 

the quantum state of the electromagnetic field, but it may depend on the boundary 
conditions, giving rise to the Casimir force O EH HH3 I13U] . 

In quantum optics, the mode operators a k are frequently represented in terms of 
the quadrature operators q k and p k |ll()j 

1 1 i 

Ofc = ^(Qk + ipk) , Qk = -j= (4 + h k) , p k = -j={a\ - a k ) . (2.10) 

The quadratures play the role of the real and the imaginary parts of the mode 
amplitudes. They satisfy the Heisenberg commutation relation (with h = 1) 

[qk,Pk>] = i$kk> ■ (2- 11 ) 

The q quadrature appears as the position and the p quadrature as the momentum 
of the electromagnetic oscillator represented in a single mode of light. This corre- 
spondence between light amplitudes and canonically conjugate quantities has found 
interesting applications in simultaneous measurements of position and momentum 
HUH [Ml UBS! and in quantum-state tomography [HU 1103 HUB 1123 HUB E22|, be- 
cause the quadratures can be measured with high precision in balanced homodyne 
detection [H 11%] . for further details see Ref. |ll()j . 



2.3 Transfer matrix 

Optical instruments act primarily on the classical wave- like properties of light. In 
the regime of far-field optics, the instrument is spatially well separated from the 
light sources and from the places where the light is detected or otherwise applied to. 
In this situation, we can decompose both the incident light and the outgoing light 
into plane waves. The reflection and transmission coefficients of the incident plane 
waves characterize the performance of the instrument. The coefficients constitute 
the transfer matrix of the dielectric structure. 
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Consider monochromatic light of frequency to propagating in a one-dimensional 
lossless dielectric. In a region where the dielectric functions e and fi do not vary, the 
mode function A k is a superposition of waves traveling to the right, exp(ikx — iut), 
and waves traveling to the left, exp(— ikx — iut), where k denotes the wavenumber 

k = ^-^ejl. (2.12) 

Consider 

A± k = ~ exp ( ^fitp) \A k ± i d x A k J , <p = j kdx. (2.13) 

In the region where the dielectric is uniform, A +k picks out the coefficient of the 
right-moving component of A, whereas A_ k gives the left-moving part. When the 
dielectric functions vary, the A± k serve to identify how the wave coefficients are 
transferred across the dielectric structure from a region of asymptotically constant 
e l and /ii on the left to a region of (possibly different) Er and /i# on the right. We 
obtain from the wave equation ()2.3|) 

/ A_ k \ _{%Z)_( 1 -exp(+2^)Ui_ fc \ 
° x \A +k - 2Z 1-expHM 1 \A +k r [ZAV 



where Z denotes the impedance [73] 



(2.15) 



In order to get reflection the impedance must vary, as we see from Eq. ()2.14jl . In 
the case of perfect impedance matching Z remains constant, and the structure is 
guaranteed to be reflectionless, a result well known from the physics of transmission 
lines [73] . To get strong reflection, with, in the extreme case, total reflection caused 
by photonic bandgaps [HIj, the dielectric structure should periodically vary at about 
twice the wave length of light, as we infer from the oscillating terms in Eq. (J2.14j) . 
We express the general solution of the differential equation (|2.14j) as 

where T(x2,Xi) denotes the transfer matrix from x\ to x%. The columns of the 
matrix T(x, x±) are required to satisfy the differential equation (|2.14|) with the initial 
condition T(xi,x\) = 1. The transfer matrix has the structure 

because, if (a, b) T solves Eq. (|2.14|) so does (6*, a*) T . Furthermore, the spatial deriva- 
tive of the determinant of T equals the spatial derivative of the impedance Z. Con- 
sequently, 

\a\ 2 - \b\ 2 = §||| , (2.18) 
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and hence 

rl = ZM t a- -6- ^ 



Z(n) V 

The transfer matrix T(+oo, — oo) characterizes the far-field performance of the op- 
tical instrument made of the dielectric structure. On the other hand, the transfer 
matrix does not directly describe how the two incident light modes interfere with 
each other to produce the outgoing modes. 

2.4 Scattering matrix 

In one spatial dimension, the directions of light propagation are fairly restricted 
- the incident light can come from the left or from the right of the dielectric 
structure. Waves incident from the left, A™, are partially reflected and partially 
transmitted, but beyond the structure the waves must propagate to the right. (We 
drop the mode index k for simplicity.) Similarly, waves coming from the right, A\ n , 
are purely outgoing towards the left. We assume monochromatic modes 

A? = ui(x) e' iujt , A 2 n = u 2 {x) e - *"* , (2.20) 

and utilize the inverse transfer matrix (j2.19j) to define u\ as the spatial component 
of a wave with the asymptotics 



Ui{X) 



\ Zl [ae ) . x^ oo _ 

L e lKRX : x — >■ +oo 



Similarly, we apply the transfer matrix (j2.17|) to define 

u 2 (x)~ A I ae -ik RX + be ik R * ; x ^ +OQ ■ (2-22) 

We normalize the A™ and A 2 n modes according to the scalar product ()2.4j) . adopting 
the procedure [HE] for normalizing Schrodinger waves in the continuous part of the 
spectrum, which also shows that the A 1 ^ and A 2 n are orthogonal to each other. We 
find 

A = 7 5 L , , 2 , M = . ? ■ |2 ■ (2.23) 

In this way we have defined the two possible incident modes for each frequency com- 
ponent of light in our effectively one-dimensional situation. Consider the outgoing 
modes. They are simply the incident modes traveling backwards, 

A° ut = u*(x) e- lujt , A° ut = u\{x) e~ iuJt , (2.24) 

forming an orthonormal set of modes as well. Since the wave equation (J2.3j) is of 
second order, each set of modes establishes a basis. Consequently, the outgoing 
modes are a superposition of the ingoing ones, 

/lout \ / Am \ 

%A=b{%), (2.25) 
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with the constant matrix B_, the scattering matrix of the beam splitter. We use the 
asymptotics (|2.21j) and (J2.22)) and the relation (|2.18j) to determine the coefficients 
of B, 

jL^) = i(VW^w -» v (2 . 26) 



a 



b* y/Z^Zl J a \ b* y/\a\ 2 -\b\ 



The scattering matrix is unitary 



Br 1 = B} . (2.27) 



Therefore, as a consequence of the mode expansion (|2.fij) . the mode operators are 
transformed in precisely the same way as the mode functions 



<;-^UU' (2 ' 28) 

To summarize this section, we may employ two alternative mode expansions ()2.6|) 
of the electromagnetic field, expansions in terms of incident or of outgoing modes, 
with the scattering matrix as mediator. The two sets of modes are adapted to two 
distinct physical situations — the incident modes refer to quantum light that enters 
the dielectric structure from outside, whereas the outgoing modes are the ones that 
leave the structure. 



3 Quantum- optical networks 

The scattering matrix completely characterizes a perfect piece of dielectric structure 
in the regime of far-field optics, describing how incident light beams are transformed 
into outgoing beams. In one spatial dimension, two monochromatic modes interfere 
to produce two emerging modes. In general, and certainly in the three-dimensional 
real world, infinitely many incident modes give rise to equally many outgoing modes. 
In experimental quantum optics, one often tries to operate with as few modes as 
possible. Much care is spent on aligning the equipment to make sure that most 
of the quantum light of interest is indeed captured in a few well-controlled modes. 
On the other hand, one can construct, in a controlled way, optical networks, also 
called multiports [188J, from the basic building blocks such as beam splitters and 
mirrors |128t 11561 1177( I192j , networks with interesting quantum properties, in 
particular in the limit when many elements are involved |176| 1178) . Furthermore, 
we could add phase-conjugating mirrors [168J or parametric amplifiers |168j to our 
catalogue of simple optical instruments, although they are experimentally less simple 
than dielectric structures. Parametric amplifiers il68j and phase-conjugating mirrors 
[168J are active devices — they require a source of energy, mostly light of a higher 
frequency. Yet these active devices share a key property with the passive instruments 
— they are linear devices in the sense that the input modes are linear transformations 
of the output modes. However, this linear transformation may involve the Hermitian 
conjugated mode operators. A phase-conjugating mirror, for example, produces the 
complex-conjugated image of the incident wave front Ak, which, in quantum optics, 
is associated with the Hermitian conjugated mode operator a\, the creation operator. 
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3.1 Linear mode transformations 



Assume that the set of mode operators at and their Hermitian conjugates, a[, de- 
scribing the incident quantum light, is turned into the operators a' k and a'l of the 
outgoing modes, by the linear transformation 




The columns (a^, a\) T and (a' k , a k ) T refer to the total set of mode operators involved 
in the transformation. We require that the operators of both the incident and the 
emerging modes are indeed proper annihilation and creation operators, subject to 
the commutation relations (|2.8jh written in matrix notation as 

\ a k J V Y a ki a k>\ l a k' a k>\ J 

with 

We substitute the mode transformation (|3.1j) and its Hermitian conjugate into the 
equivalent relation for a' k and a k , and find 

SGS} = G. (3.4) 

Such transformations are called quasi-unitary [IS]. In the classical mechanics of a 
many-particle system they are called linear canonical transformations, because 
they preserve the canonical Poisson brackets (HHl- Equation (|3.4|) implies that G_S}Q_ 
is the inverse of S_. Therefore, S_ is a square matrix. We get for the determinant 

I det 5; | 2 = 1 . (3.5) 

Consequently, the Jacobian of the mode transformation has unity modulus — the 
phase-space volume is conserved, as we would expect from canonical transforma- 
tions according to Liouville's theorem j^Hl- Quantum mechanics requires that the 
number of input modes is exactly the same as the number of output modes. Modes 
that are "empty" contribute nevertheless to the quantum properties of the device. 
They are not really empty, they are just in the vacuum state. The vacuum noise 
behind a mirror matters and so do the vacuum fluctuations of modes prior to 

amplification. 

When the optical instrument transforms annihilation operators into annihilation 
operators without involving their Hermitian conjugates, as it is the case for the one- 
dimensional dielectric structures analyzed in Sec. 2 or for passive optical multiports 
in general [US ESI ESI E3 EFJ , we get 

- = ( f B* ) ' -- = 1 • (3 ' 6) 
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We see that the unitarity ()2.27|) of the beam-splitter matrix B_ is not a coincidence - 
it follows from the conservation of the commutation relation (J2.8j) in light scattering. 
As a consequence of the unitarity of B_, the total number of photons is conserved, 

k k 

Passive optical instruments conserve the total energy. Active devices such as phase- 
conjugating mirrors |168j or parametric amplifiers [168J combine annihilation and 
creation operators. Consequently, the total number of photons is not conserved, in 
general, indicating that active devices rely on external energy sources. 



3.2 Quantum-state transformations 

The linear transformation (j3.1|) of mode operators describes the transformation of 
the incident quantum light into the emerging quanta in a peculiar manner. In the 
mode transformations, the quantum state of light is invariant, but its relation to 
physical observables changes, similar to the Heisenberg picture of quantum mechan- 
ics. A given quantum superposition of photons, frozen in space and time, is seen 
first as constituting the incident modes and then as leaving in the outgoing modes. 
The Schrodinger picture gives perhaps a more natural approach to understanding 
the quantum effects of optical instruments on light. Here the instrument changes 
the state of light, whereas the mode operators remain invariant. In other words, in 
the Schrodinger picture the incident and the emerging modes are the same and the 
optical instrument operates like a black box on the quantum state of light. This 
picture is especially suitable for analyzing laboratory situations where a few well- 
controlled modes enter the instrument and leave it in other equally-well-controlled 
modes, describing how the quantum state is transformed. The Heisenberg and the 
Schrodinger picture ought to agree on their quantum-mechanical predictions, on 
expectation values, and this is how we deduce the quantum-state transformations 
from the linear mode transformations (|3.1|) 118j . 

We describe the quantum state of light in terms of the density operator (also 
called density matrix) j!2l EH HUH IHUj • We require for any physical observable of 
light, for any function of the mode operators, that the expectation values in the two 
pictures agree, 

tr{p/(a' fc , a'l)} = tr{p'f(a k , a\)} , (3.8) 

where p' denotes the density operator of the outgoing light. If we manage to deter- 
mine a unitary evolution operator B with the property 




we get 

p=B ] pB. (3.10) 
Consider the logarithm of S_ defined as a matrix In S_ for which 

expQnS) = V = S. (3.11) 

n=0 
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The quasi-unitarity ()3.4j) of the S_ matrix implies that — In 5 = ln(S_ 1 ) = \n(GS}G). 
Furthermore, since G 2 = 1, we obtain from the definition (J3.ll)) of the matrix 
logarithm that \n(GS}G) = G (ln5) f G. Consequently, 

H = -iG\nS (3.12) 

is a Hermitian matrix. We construct the operators |118j 

B = exp(iH), H = ±(al,a k )H( a £). (3.13) 

Since H is Hermitian, B is unitary. We prove that B does indeed act as an evolution 
operator with the property (|3.9jl . Consider the power i? 7 ' for real 77. First we show 
that the differential equation in 77 of B v akB^ v is the same as the differential equation 
for the transformed mode operators with matrix S_ v . From the commutation relations 
()2.8J1 in matrix form ()3.2j) follows 

= iexp(ir]H) GH_( a \ j exp(-ir)H) 
= iGHB n (J*) (3.14) 
which indeed agrees with the differential equation 

d v S*( ^)=iGHS!>( ^) (3.15) 



for if given in terms ()3.12|) of the matrix logarithm. Since at 77 = the mode 
operators are not transformed, the initial B v akB^ v trivially agrees with the effect of 
S_° = 1. Therefore, B v (a,}., a\.) T B^' n gives S n (ai t ,a1) T all the way up to r\ = 1, thus 
proving the relation ([3.9)1 . 

The H operator plays the role of the effective Hamiltonian for the quantum- 
optical network |175| I177j , generating the linear mode transformation ()3.1)) in the 
Heisenberg picture. This Hamiltonian depends on the logarithm of the transfor- 
mation matrix 5, which is a multivalued function with infinitely many branches. 
So there are many equivalent ways to design an optical network with a particular 
input-output relation and there are also many ways to assemble it from the basic 
building blocks [156. 1175) 1177] . from beam splitters and parametric amplifiers. 

3.3 Examples 

The two prime examples of simple optical instruments are the beam splitter and 
the parametric amplifier. We describe the beam splitter by the unitary scattering 
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matrix B_. For simplicity, let us assume that B_ is a real rotation matrix instead of 
the general 2x2 unitary matrix, 



B 



cos if sin if 
— sin if cos (f 



(3.16) 



with the rotation angle if. For example, the matrix ()3.16j) may describe a polarizing 
beam splitter where an incident light beam is separated into two linear polarizations 
with angles if and if + n. We show in Sec. 4 that the simple rotation matrix 
(|3.16|) contains the essence of all two-mode beam splitters. Here we derive the 
Hamiltonian for the linear mode transformation (|3.1|) described by the beam-splitter 
matrix (|3.16|) . Consider the matrix 



(3.17) 








1 





o\ 




-1 






















1 










-1 





Since V 



•1 we get 



k=0 



2 k 



2fc+l 



/c=0 



(2k + l)\ 



-i) k i 



= (cos^)l + (simp) I. (3.18) 

Consequently, ifl_ is the logarithm of the S_ matrix with S_ given by Eq. (J3.6j) and 
the beam-splitter matrix ()3.16|) . Therefore, the effective Hamiltonian ()3.13|) of the 
beam splitter is 

H — iif (a\a 2 — o^ij • (3.19) 

The Hamiltonian indicates that photons from mode 1 are annihilated and converted 
into photons of mode 2, and vice versa, which is just what we expect from a beam 
splitter. As we know, the total number of photons is conserved for passive optical 
instruments. 

Let us turn to active devices that, by definition, mix annihilation and creation 
operators. A simple example of an active device is characterized by the S_ matrix 



S 



( cosh ( 



V sinh£ 





cosh ( 
sinh£ 






sinh ( 
cosh£ 




sinh£ \ 



cosh£ J 



(3.20) 



with the real parameter (. One easily verifies that the matrix ()3.2()j) indeed satisfies 
the quasi-unitarity relation ()3.4j) and thus qualifies for the S_ matrix of an optical 
instrument. In fact, the matrix describes a parametric amplifier |168j or a phase- 
conjugating mirror |168j . see Sec. 6. We derive the effective Hamiltonian. Consider 
the matrix 



E 



(0 








1 \ 








1 








1 








V 1 











(3.21) 
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Since Ej = 1 we get 

exp«£) = (cosh C)l + (sinh ()E = S . (3.22) 
Consequently, the effective Hamiltonian ()3.13|) is 

H = iC, (^a[al — &iQ>2 J ■ (3.23) 

The Hamiltonian indicates that two photons are simultaneously created or annihi- 
lated. The pump process of the amplifier, accounted for in the parameter (, must 
provide the energy source of the photon-pair production or the reservoir for annihila- 
tion. We could represent ( as '"ft/2 where t denotes the amplification time and 7 the 
differential gain that depends on the performance of the pump. In our simple model 
(|3.2(J|) the pump is assumed to be classical and to remain essentially unchanged, 
giving rise to a constant rate 7 during the amplification. 



3.4 Wigner function 

Quasiprobability distributions [IUJ 11101 H61| are frequently used in quantum optics 
to draw intuitive pictures of the quantum fluctuations of light, for computational 
advantages and to give a precise meaning to the notion of non-classical light |llUj . 
Quasiprobability distributions are functions of the classical quadratures q and p 
that behave in many ways like classical probability densities. However, since the 
quadrature operators q and p do not commute, since q and p cannot be measured 
simultaneously and precisely, the quasiprobability distributions must not represent 
perfect phase-space densities. For example, they may appear to describe negative 
probabilities or they may become mathematically ill behaved |llUj . In quantum 
optics, the most prominent quasiprobability distributions are the P function, the 
Q function and the Wigner function, see for example Ref. |110j . The P function 
allows us to express, by the optical equivalence theorem [oT)lllimll71| . any quantum 
state as a quasi-ensemble of coherent states (of classical light waves) |llUH123lll27j . 
If the P function is non-negative and well-behaved the light is said to be classical. 
Otherwise the light is non-classical in the sense that it cannot be understood as 
partially coherent classical light. The Q function is proportional to the expectation 
value of the density matrix in coherent states [110j . The Q function appears as 
the genuine probability distribution in simultaneous measurements of position and 
momentum quadratures |llUU187lll88| . However, since q andp cannot be measured 
both simultaneously and precisely, the Q function contains some extra quantum 
noise that is difficult to remove from the true quantum state by deconvolutions 

[IHl EDI- 

The Wigner function |HUt 11341 11611 1173| I193j is probably best suited to describe 
the quantum effects of simple optical instruments. The Wigner function W(q,p) of 
a single mode of light is the inverse Fourier transform of the characteristic function 

w(u,v) mm 

W(q,p) = - — — / / W(u, v) exp(iuq + ivp) du dv , 
(2vr) 2 y_ 00 

W{u,v) = tr{pexp(— iuq — ivp)} . (3.24) 
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We obtain |llUj in terms of the q quadrature eigenstates (position eigenstates) 



-y r+oo 

W(q,p) = — / exp(ipx)(q — x/2 | p \ q + x/2) dx . (3.25) 

The Wigner function is real and normalized to unity for any proper density operator 
110J. Quantum expectation values can be computed via the overlap formula |llUj 

/+oo r+oo 
/ W 1 (q,p)W 2 (q,p)dqdp, (3.26) 
-oo J — oo 

where W\ and W 2 are the Wigner transforms with p in formula ()3.25j) replaced by 
F\ and F 2 , respectively. The Wigner function gives a faithful and frequently quite 
intuitive image of the quantum state of a single mode of light. Optical homodyne 
tomography has been applied to reconstruct the Wigner function from homodyne 
measurements [HU I110| 11261 H7U| I192j . The marginal distributions of the Wigner 
function agree with the correct quadrature histograms with respect to an arbitrary 
phase shift. This tomographic principle underlies optical homodyne tomography 
and it also defines the Wigner function uniquely [2U IllOj . The Wigner function 
represents a fairly good compromise between the abstract density operator of quan- 
tum mechanics and the phase-space density of classical statistical mechanics, but 
the Wigner function may exhibit negative "probabilities" in small phase-space re- 
gions |ll()j . Such features are quite subtle and have been observed only recently in 
quantum light |231IT2B| . 

Here we use the Wigner function to describe the quantum effects generated by 
simple optical instruments that are subject to the linear mode transformations (|3.1|) . 
We extend the definition (|3.24|) of the Wigner function to a multitude of light modes 
characterized by the classical amplitudes 

a k = -j={Q.k + Wk) ■ (3.27) 

We represent uq k + vp k in the definition ()3.24|) of the characteristic function as 
Pla k + (3 k a\ with (3 k = (u k + iv k )/\/2. We see that 



tr jp' exp f-i Y^(P* k a k + ^ fe oJ.)j 

tr jp exp ^-i Y^(f3* k a' k + (3 k a k f )j 
WW, ft) (3.28) 



^ ^ - ST 1 ( ^ ) . (3.29) 



with 

To obtain the Wigner function ()3.24|) of the emerging multi-mode light we represent 
uq k + vp k as Pla k + (3 k a* k in the inverse Fourier transformation of the characteristic 
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function and we use fa and (51 as the integration variables. Then we perform a 
variable transformation from fa and fa to fa and fa k *. The Jacobian (|H.5jl of this 
transformation has unity modulus, and we get the result 

W'(a k , al) = W(a' k , a';) , ( < ) = ST 1 ( ^ ) ■ (3-30) 

Simple optical instruments transform the Wigner function W of the incident quan- 
tum light as if W were a classical probability distribution of the mode amplitudes. 
This property uniquely distinguishes the Wigner function for general quasi- 
unitary transformations involving Hermitian conjugated mode operators, i.e. for 
active optical instruments |lU4j . Passive instruments such as optical multiports 
156 Ul75| I177| transform also the P function and the Q function like classical prob- 
ability distributions |lU2j . 

4 Beam splitter 

The archetype of passive optical instruments is the beam splitter, usually an innocent- 
looking cube of glass in laboratory experiments, see Fig.^ Light incident at the front 
of the cube is split into two beams, and so is light incident at the back. Both modes 
may interfere. In Sec. 2 we studied the theoretically simplest example of a beam split- 
ter, a one-dimensional dielectric structure. Polarizers, where the two polarization 
modes of light are mixed, are also essentially beam splitters, and so are simple passive 
interferometers. The theory of this passive four-port device has been developed in 

Refs. m ehi eh IZ2J cznj izzi izhj E3i [fn2i H2H hsbi nisi mni unoi nnsi nHsi nnsi - Here 

we follow mostly Refs. |H |ll02j . More complicated optical multiports jl56|IT75]ll77j . 
where a multitude of beams interfere to produce the same number of outgoing modes, 
can be constructed from beam splitters and mirrors |156j . But already the simple 
beam splitter, combined with good photodetectors and single-photon sources, is 
quite capable of demonstrating some fundamental aspects of the wave-particle du- 
ality of light. 

4.1 Matrix structure 

The beam splitter is completely characterized by a unitary 2x2 matrix B_ that 
describes how the device transforms the incident modes into the outgoing modes in 
the Heisenberg picture, 

The beam-splitter matrix is unitary, in order to preserve the Bose commutation 
relations between the mode operators. Explicitly, the matrix elements must obey 

I Bu \ 2 + | B12 | 2 = 1 , | B 2 ± | 2 + | B 2 2 | 2 = 1 , Bl x B 2 i + Bl 2 B 2 2 = . (4.2) 
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a 2 

Figure 1: Beamsplitter. The picture above shows a polarizing beam splitter. (Cour- 
tesy of Oliver Glockl and Natasha Korolkova.) The picture below schematically 
illustrates the quantum theory of the beam splitter. Two incident light modes, rep- 
resented by the Bose mode operators Si and 02 interfere to produce two emerging 
modes with operators d[ and a' 2 . Even if only one incident beam is split the vacuum 
noise of the second incident mode behind the mirror plays an important role in the 
quantum optics of the beam splitter. 
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The general solution of these equations is 



iA/2 f cos(6/2) e*^*)/ 2 sin(6/2) e^"*)/ 2 \ 
-~ 6 V -^(e^e^-*)/ 2 cos(e/2)e-^ +,I> )/ 2 y< l4 ^ j 

with the real parameters $, 0, \1/ and A. The one-dimensional dielectric structure 
with matrix ()2.26|) represents the special case where $ = — We express the 
general beam-splitter matrix ()4.3|) as the product 

_ piA/2 / e**/ 2 \ / cos(0/2) sin(0/2) \ / e^ 2 \ 
-~ 6 V e_< * /2 / V - sin(0/2) cos(0/2) J \ e^' 2 ) ' 

The beam splitter acts in four steps. The incident modes gain a relative phase of $, 
the modes are optically mixed with the mixing angle 0/2, and the outgoing modes 
attain the relative phase \1/ and the overall phase A/2. We could incorporate the 
phases into the definitions of the incident and the outgoing modes. The rotation 
matrix would remain as the key feature of the beam splitter. The reflectivity g is 
characterized by — sin(0/2) (the sign is unimportant though) while the transmis- 
sivity t is given by cos(0/2), which implies r 2 + g 2 = 1. The beam splitter has 
been assumed to be perfectly lossless — if a photon is not transmitted it must be 
reflected. We show in Sec. 5 how, in principle, absorption can be included and that 
the beam splitter itself serves as a convenient model of an absorber. 



4.2 Quantum Stokes parameters 

The classical polarization of a light beam is usually described using the Stokes pa- 
rameters (27J- Given the complex amplitudes a\ and a 2 of the two polarization 
modes, the three Stokes parameters are proportional to the corresponding expecta- 
tion values of the Pauli matrices 

a * = (l I ) ' °v = ( i 7 ) ' ° z = ( I -1 ) ■ (4 ' 5) 

The Stokes parameters lie on a sphere, the Poincare sphere [27], also called the Bloch 
sphere in quantum mechanics |127j , see Fig. El In quantum optics, we describe the 
polarization, the spinor part of the angular momentum of light, in terms of the 
quantum Stokes parameters 
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Figure 2: Poincare sphere. The classical polarization of light is characterized by 
the Stokes parameters that lie on the Poincare sphere. In quantum optics, the 
Jordan- Schwinger operators (J4.6|) play the role of the quantum Stokes parameters 
and their quantum statistics characterizes the polarization state. The parameter L t , 
the conserved total photon number, corresponds to the radius of the Poincare sphere. 
The polarizer of Fig. ^rotates the quantum Stokes parameters. In general, we can 
characterize the two incident modes of any beam splitter by their Jordan-Schwinger 
operators as well, and the beam splitter performs a rotation on the generalized 
Poincare sphere. 
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that obey the commutation relations of angular-momentum operators 



[L X , Ly] 



iL, 



[Ly, L 2 



iL n 



[L z , L x ] 



iL 



y ■ 



(4.7) 



Equation (J4.6)) is called the Jordan- Schwinger representation |HHlll66j of the angular 
momentum in terms of two Bose operators. The representation serves as the starting 
point for the quantum theory of polarized or partially polarized light jHHl HOlj • The 
operator L t commutes with all others and serves to represent the squared total 
angular momentum, 

L 2 x + L 2 y + L 2 z = L t (L t + l). (4.8) 

The commutation relations (J4.7)) give rise to uncertainty relations between the quan- 
tum Stokes parameters. Polarization squeezing [HH] occurs when the statistical fluc- 
tuations of one Stokes parameter are below the minimum-uncertainty limit. This 
quantum-noise reduction of the polarization of light has been applied to observe 
macroscopic spin-squeezing effects in atomic vapors [BB1 E3j • 

The Jordan- Schwinger representation 1)4.6)1 serves not only to characterize the 
polarization of quantum light, the representation provides also the theoretical tools 
to describe the effect of polarizers or of any beam splitter in general. We obtain 
from the commutation relations (J4.7)) 



exp(— z$L- 




exp(— iflL x ) 



L, 



exp(z$L 2 



L x 

exp(-iQLy) | L y | exp(i®L y ) 



exp(iQL x ) 




(4.9) 



as one easily verifies by differentiation with respect to the parameters. Therefore, the 
exponential Jordan- Schwinger operators describe rotations on the Poincare sphere, 
generated by polarizers. We note that the angles $, 0, \I/ in the complex matrix 
representation (J4.3)) are the Euler angles of an arbitrary rotation in three-dimensional 
space j^Hl- We obtain for the mode operators 
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in agreement with our previous result 1)3.19)1 for the effective Hamiltonian of the 
real beam splitter (jM.lfij) . Consequently, the beam splitter performs rotations in the 
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three-dimensional space spanned by the quantum Stokes parameters ()4.6|) . Such 
rotations are independent on the overall phase A/2 in the factorization (J4.4j) . For 
rotations on the quantum Poincare sphere we can thus restrict the beam-splitter 
transformations to SU(2) matrices jlHI with 

detB = l. (4.11) 

Finally, we arrive at the general evolution operator for the beam splitter with matrix 

B = exp(-z$L z ) exp(-i9Lj exp(-WL z ) exp(-iAL t ) . (4.12) 
4.3 Wave-particle dualism 

Now we possess the theoretical tools to predict what happens when two beams of 
quantum light interfere at a beam splitter. In classical optics, the light beams are 
characterized by their spatial shapes, by their normalized spatial mode functions U\ 
and u 2 , and by their amplitudes a x and a 2 . The outgoing modes have the amplitudes 




(4.13) 



The amplitudes of the incident light modes may statistically fluctuate if the beams 
are not perfectly coherent [23 I127j . In this case the outgoing modes fluctuate 
accordingly, because the individual amplitudes are transferred according to the re- 
lation (|4.13|) . In quantum optics, the observables such as the amplitudes di and 
a 2 or the photon numbers a\ai and a\a2 may statistically fluctuate in repeated ex- 
periments, even if the light has always been prepared in identical pure states |llUj . 
Such quantum fluctuations tend to be quite subtle and hard to discriminate from 
classical noise in experiments. The quantum-noise properties distinguish the vari- 
ous quantum states of light |ll()j . The coherent states |ll()j resemble classical light 
beams with well-defined amplitudes, coherent light. A coherent state of light is 
characterized by the Wigner function |110j 

W(q,p) = - exp (-{q - q ) 2 - {p - p ) 2 ) = - exp (-2| a - a \ 2 ) , (4.14) 

71 71 

where a = (q + ip )/\/2 denotes the classical amplitude of the light beam and 
a abbreviates (q + ip) / y/2. The Wigner function ()4.14|) describes the quantum- 
statistical fluctuations of the amplitude components q and p, the quadratures |ll()j . 
The vacuum state belongs to the class of coherent states as well [110J — the vacuum 
is the coherent state with zero average amplitude. Yet the field amplitudes of the 
vacuum state still fluctuate |f 10j . We see from the Wigner function (j4~T4j) that the 
amplitudes of coherent states fluctuate precisely like the quantum vacuum around 
their average values a,Q. The coherent states are the most classical-like states of 
light, corresponding to waves as perfect as quantum mechanics allows. Note that 
the number of photons fluctuates in a coherent state, because precision in amplitude 
and precision in particle number are mutually exclusive. The photons in a coherent 
state of amplitude qq are as randomly distributed as raisins in a cake with | ao| 2 
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apiece |llUj . Technically |llUj . the photons follow a Poisson distribution around the 
average | a 1 2 . 

Suppose that the two light beams incident on the beam splitter are in coher- 
ent states with the amplitudes a\ and a 2 , corresponding to the two-mode Wigner 
function ^ 

W(ai, a 2 ) = — r exp (—21 oi\ — a x \ 2 — 2| a 2 — a-2 1 2 ) • (4-15) 

To predict the quantum state of the outgoing modes, we apply the transformation 
rule ()3.30|) to the Wigner function (j4.15j) and utilize the unitarity ()2.27|) of the 
beam-splitter matrix. We obtain 

W'(a 1 ,a 2 ) = At exp (-2\ ai - a[\ 2 - 2| a 2 - a' \ 2 ) . (4.16) 

The outgoing modes are in coherent states with the amplitudes classically trans- 
formed (|4.13|) . The modes are completely uncorrelated, because the Wigner function 
factorizes. The coherent states thus interfere just like classical waves, even down to 
the finest details of their quantum-statistical properties. This result uniquely dis- 
tinguishes coherent states [I] and it can be extended to any passive optical network 
|156t I175| 1177] . Historically, the interference property of coherent states has been 
deduced from a microscopic model of the beam splitter 38] and has served as the 
starting point for the quantum theory of such optical instruments |143j . 

Now, suppose that one incident beam carries precisely n photons and that no 
light impinges on the back of the semi-transparent mirror. The light beam with 
exactly n photons is in the Fock state jllOj 

|n) = ^a tn |0), (4.17) 



and the other incident mode is in the vacuum state | ) . Fock states are the eigen- 
states of the photon-number operator a) a and hence they correspond to light with 
a perfectly well-defined number of photons |llUj . We calculate the quantum state 
of the outgoing modes 

£t| n ) x I ) 2 = -L £t a t»| 0>i I 0> 2 = A= fraPB I > a I > 2 - (4.18) 
yn\ vn! 



Here we have used the fact that the beam splitter transforms the incident vacuum 
into the outgoing vacuum, ex nihilo nihil, as we easily see from the Wigner function 
(I4.16|) with zero initial amplitudes a\ and a 2 . Since 

we obtain according to the Binomial theorem and the definition (I4.17J1 of the Fock 
states 

B^\n) 1 \0) 2 = -= (B n a\ + B 21 aU \0h\0), 




B^B^lkhln-kh, (4.20) 
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The beam splitter does not split the incident photons, of course, but rather the 
semi-transparent mirror statistically distributes the photons into the reflected and 
the transmitted beam. Suppose we count the photons in each emerging mode |3*2*] . 
Each individual run of the experiment [32] is unpredictable, but averaged over a 
large statistical ensemble we get the joint photon-number distribution 

n \ Jlk /i J2\n-k 



p(ni,n 2 ) = ( I ) T 2k (l-T 2 ) n - k 5 nuk 6 n2 , n . k} (4.22) 

where r denotes the transmissivity \Bu\. The Binomial distribution ()4.22j) describes 
a random decision process where n distinguishable objects are distributed to two 
channels, to the first channel with probability r 2 per object and to the second one 
with probability 1 — r 2 , accordingly, because the beam splitter is assumed to be 
perfectly lossless. Each photon is statistically independent, and so the probability 
for k individual photons to arrive in the first channel and n — k photons in the second 
one is the product r 2k (1 — r 2 ) n_fc . We multiply this value by the Binomial coefficient, 
which describes the number of possibilities to distribute any k of the n photons to 
the first channel and the rest to the second one, because we cannot discriminate 
between individual photons in photon-counting experiments. Nevertheless, photons 
behave in beam-splitting experiments as if they were in-principle distinguishable, in 
contrast to the common statement that photons are fundamentally indistinguishable 
particles, which illustrates some of the conceptional subtleties of the photon. 

There is another twist in the physics of photons and the beam splitter. Suppose 
you let two beams of light with equal intensities interfere at a perfect 50:50 beam 
splitter characterized by the real matrix 

Consider two coherent states with equal complex amplitudes CLq — <X\ — CL2- We 
obtain from the transformation rule (j4.13|) that the first outgoing mode is in the 
coherent state with complex amplitude \/2ao, whereas the second mode is in the 
vacuum state. The two incident light beams interfere constructively in the first 
outgoing mode and destructively in the second one. If the two incident coherent 
states have equal amplitudes but opposite phases, a\ = Oq, 02 = — a Q , they interfere 
the other way round. Now, suppose you let one single photon interfere with another 
single photon. We calculate the quantum state of the outgoing modes, 

£ + |i)i|i) a = ^(4-4) (4 + 4) |o>i|o) 2 

= -L(|2) 1 |0) 2 -|0) 1 |2) 2 ). (4.24) 

The photons interfere constructively or destructively. Complete destructive inter- 
ference implies that the affected outgoing mode is in the vacuum state, whereas the 
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other mode must carry exactly two photons, because the total number of photons is 
conserved. Destructive or constructive interference depends on the relative phase of 
the incident photons. Fock states with precisely defined photon number are the most 
extreme particle-like states of light, and hence they do not carry any wave-like phase 
information. Faced with this dilemma, the beam splitter distributes the two pho- 
tons in either way with l/y/2 probability amplitude, i.e. with 50% probability after 
detection, a truly Solomonic solution. Experimentally |7T] . the coincident counts of 
photons reach a well-pronounced minimum when the spatial-temporal modes of the 
incident photons overlap at the beam splitter. The choice which one of the outgoing 
mode carries the two photons becomes only apparent when the light is detected. The 
same is true for the beam splitting of n photons discussed previously. The beam 
splitter itself is a deterministic device. The probabilistic outcome of the photocount- 
ing remains undecided until the measurement is made. The quantum state (J4.24j) 
of the outgoing modes is strongly correlated, and so is the state ()4.20|) of the split 
Fock state, in contrast to the interference ()4.16|) of coherent states. Moreover, the 
decisive measurement devices may be located a long distance apart from each other. 
Both photon-interference and photon-splitting experiments are suitable |172^ IMj to 
test the non- locality of quantum mechanics fK)\ I5T] . 

5 Absorber and amplifier 

The optical instruments studied so far are completely reversible devices. For exam- 
ple, when a Fock state, carrying a precise number of photons, is split at a beam 
splitter we could, in principle, send the outgoing beams back to restore the initial 
Fock state. Mathematically argued, the instruments (|3.1|) are reversible, because 
the quasi-unitary S_ matrix has the inverse GS}G, according to Eq. (|3.4|) . So far, 
we have excluded irreversible processes such as the absorption or the amplification 
of light. However, all quantum processes are fundamentally reversible, as long as 
no measurements are made or could be made in principle (whatever measurement 
processes are), and as long as we keep track of the quantum systems involved. 

Consider for example an absorber, a piece of grey material. Some of the incident 
light is destined for absorption and some part is transmitted, with reduced intensity 
though. The absorbed component is transferred to the material and ebbs away 
in many small material excitations as heat. Normally we are simply not able to 
keep track of the material details and so the absorbed quanta are lost. (Exceptions 
are very simple atomic systems, two-level systems for example, where absorption 
is reversible [7^ 11421 H69| .) The first part of the absorption process resembles 
beam splitting. We could model the second, the irreversible part by discarding the 
quantum information carried in one of the outgoing beams, i.e. by averaging over 
the unobserved component of the total quantum system. Modeling absorption and 
detection losses by fictitious beam splitters has been a successful idea in quantum 
optics jSHl EDI IHEl I103| 11451 1194j . This theoretical trick is known elsewhere as the 
thermo-field technique I 1 80: . In fact, many absorbers are first of all scatterers, but it 
is quite remarkable that we can sum up the multitude of scattered light modes in just 
one outgoing mode of a fictitious beam splitter. Moreover, an absorber may emit 
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thermal radiation according to its temperature, and so we should include emission 
as well as absorption in our model. When the stimulated emission dominates the 
device acts as an amplifier. To understand how to model such irreversible processes 
requires some theory [33J E21 EH GUI H20| . 

5.1 Lindblad's theorem 

Lindblad [120J determined the most general structure of the dynamic equation for 
the density operator p, assuming only that the evolving p represents indeed an 
ensemble of pure quantum states | ip a ) occuring with probabilities p a . A reversible 
quantum process would only change the pure states | ij) a ) while an irreversible process 
may affect both the | ip a ) and their probabilities p a . Lindblad's master equation, the 
quantum version of the Boltzmann equation, reads [SSI 132 EHJ E3 H20j 

d t p = l - [p, Ho] - X> (UUP ~ 2LipL\ + pL\L^ . (5.1) 

Here Hq denotes the Hamiltonian of the reversible part of the dynamics, while the 7/ 
quantify the rates of the irreversible processes described by the Lindblad operators 
Li. The parameter t may describe a fictitious time that, for example, corresponds to 
the penetration depth of an absorbing material or to the length of a laser amplifier. 
We define the effective Hamiltonian 

H eS = H -ihJ27iL 1 iL l , (5.2) 
1 

a non-Hermitian operator, and write the master equation (j5.1j) as 

8 t p = % - [pH\ n - H cS p) + 2 liUpl] . (5.3) 

1 

The effective Hamiltonian alone would reduce the total quantum probability trp. 
The component i(pHl s —H cS p)/h of the master equation (|5.3|) describes the coherent 
part of the irreversible process, for example damping or amplification, while the 
terms 2^iLipL[ characterize the effect of quantum jumps |3*3*1 142j . fluctuations that 
restore the nature of the density operator. In this respect, Lindblad's theorem [120J 
formulates the quantum version of the fluctuation-dissipation theorem [99J. 

The Lindblad operators L t describe the specific physical effects of the irreversible 
processes involved in the dynamics (|5.1|) . the quantum transitions caused. In short, 
the Li are the transition operators. For specific irreversible processes we can fre- 
quently use our intuition to infer the relevant transition operators. We may guess 
that the absorption of light corresponds to the effect of the annihilation operator d, 
while the light emission is represented by the creation operator a*. An absorber or 
amplifier is thus modeled by the Lindblad operators 

Li = a , L 2 = a) . (5-4) 

For simplicity, we ignore the Hamiltonian H of the single light mode that would only 
generate a time-dependent phase shift of light. We translate the master equation 
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()5.1|) with the Lindblad operators (|5.4j) into the evolution equation of the Wigner 
function, a Fokker-Planck equation |4*2*ll5b|lo"7|ll59| . For this, we express the master 
equation in terms of the quadratures q and p with a — (q + ip)J V2 and calculate 
the Wigner transforms (|3.25|) of the operators involved. We use the correspondence 
rules jSHl 

qF^(q+ l -d)jW F , Fq^(q- l -dAW F} 

pF^[p-\d)jW F , Fp^ (p+'-d^Wp, (5.5) 

between the operators and their Wigner transforms, and arrive at the Fokker-Planck 
equation 

d t W = (71 - 72)(d q (qW) + d p ( P W)) + ?±±-^(d 2 W + d 2 p W) . (5.6) 

The first term, with prefactor 71 — 72, describes the drift of the quasiprobabilities, 
to zero if the absorption dominates and to infinity if the emission is stronger. The 
second term, with rate (71 + 72V2, describes the diffusion of the quasiprobabilities 
due to quantum noise. 



5.2 Absorber 



Suppose that the absorption rate 71 outweighs the emission rate 72. In this case, 
the Fokker-Planck equation ()5.6|) describes the net effect of an absorber. We find 
the stationary solution, normalized to unity, 
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The Wigner function W^h corresponds to the thermal state |llUj 
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with average photon number N and temperature T, according to Planck's formula 
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(5.9) 



The stationary solution indicates that the absorber consists of a thermal reservoir 
with temperature T, a reservoir that absorbs light, but that also emits thermal 
radiation. Now, consider the general solution of the Fokker-Planck equation (|5.fij) . 
One verifies easily that the solution is [102 j 
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(5.10) 
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where the a's abbreviate the (q + ip)/\/2 amplitudes, with 



a 'i \ _ ( \fn -V 1 ~v \ ( ot 



a 2 J V V 1 - V Vv J V a ° 



(5.11) 



and 

77 = exp[-2( 7l - 72 )(t-to)] > *>*o- (5-12) 

Any initial Wigner function W(a,t ) is exponentially attenuated and eventually 
approaches the thermal state W t h in the limit t — > +oo. The result (|5.1U|) with 
the relation (|5.11|) proves that partial absorption corresponds to a simple beam 
splitter model. The incident light appears to be split into the transmitted and the 
absorbed component. The rj parameter (|5.12|) describes the transmission probability 
of a single photon and 1 — rj is the probability of absorption. The second mode of 
the fictitious beam splitter plays a double role. The mode represents the reservoir 
into which the absorbed light disappears and over which we average in the solution 
f)5. lOj) . Additionally, the initial state of the mode describes the fluctuations of the 
thermal reservoir that contaminate the transmitted light. In the Heisenberg picture 
we represent the mode operator of the partially absorbed light as 

a(t) = a(t )^/rj + a y/l - V . (5.13) 

The fluctuation mode is essential in order to preserve the Bose commutation relation 
of a(t), even at zero temperature. Usually, for light in the optical range of the 
spectrum, room temperature and zero temperature makes little difference. In this 
case, the vacuum fluctuations of the reservoir affect the partially absorbed light. If 
the light has initially been in a coherent state with amplitude a the transmitted 
light remains in a coherent state with the reduced amplitude ^a , despite the 
vacuum fluctuations, because, as we know, the beam splitter transforms coherent 
states (the initial state and the vacuum mode) into disentangled coherent states. On 
the other hand, other quantum states approach coherent states during the absorption 
process. The absorber purifies light with excess amplitude noise, but the absorber 
also destroys fragile non-classical states with unusual quantum properties, such as 
Schrodinger-cat states |3*7 | 1102] . 



5.3 Amplifier 

Suppose that the emission rate 7 2 outweighs the absorption rate 7 i in the irreversible 
process (|5.1|) with the Lindblad operators (|5.4|) . In this case the Fokker-Planck 
equation ()5.6|) has the general solution ()5.10|) with [104J 



'i 

j* 



^ I VV -VV- I W " \ (5.14) 



a 2 ) V -y/V 1 y/V J \ a < 



'0 



and with the rj parameter (|5.12j) larger than unity. The amplitude a of the initial 
Wigner function grows with ^/rj and exponentially in t, which indicates that the 
process describes a linear amplifier with gain rj = exp[2( 7 2 — 7 i)(t — to)]- The ampli- 
fication is accompanied by amplification noise, usually spontaneous-emission noise 
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in laser amplifiers, summed up in the thermal Wigner function ()5.7|) that becomes 
interwoven with the initial quantum state. The noise temperature (|5.9|) character- 
izes the quality of the amplifier. Amplification noise is stronger than absorption 
noise in the sense that amplified coherent states do not remain coherent states, even 
at zero noise temperature. Therefore, amplification does not simply reverse atten- 
uation. The growth of the signal mode combined with the inevitable amplification 
noise appears in the Heisenberg picture as 

a(t) = a(t )</fj + a\^r] - 1 . (5.15) 

As in the case of the absorber, the fluctuation mode preserves the Bose commutation 
relation of the amplified light, but the amplification noise always creates additional 
quanta, indicated by the creation operator, in contrast to the absorption noise. 

Suppose that the amplifier attempts to balance the effect of attenuation, i.e. 
7i = 72 = 7 in the process (|5.1j) with the Lindblad operators ()5.4j) . Imagine, 
for example, that Eve, the eavesdropper, tries to tap quantum information by beam 
splitting while covering up her tracks by amplification. In the case when the emission 
rate is equal to the absorption rate the Fokker-Planck equation ()5.6|) reduces to the 
pure diffusion of the Wigner function, without drift, as designed. We calculate the 
quantum-statistical purity [HU IllDj of the signal state using the overlap formula 

(USED 

/+oo ^+oo 
/ W 2 dqdp, (5.16) 
-oo J — oo 

and obtain from the Fokker-Planck equation ()5.fij) by partial integration 

/+oo ;>+oo 
/ ((dgW) 2 + (d p W) 2 )dqdp < 0. (5.17) 
-oo J —oo 

The purity tr{p 2 } monotonously decreases until the Wigner function has been com- 
pletely leveled by diffusion, containing no information anymore. Eavesdropping 
spoils the purity of the quantum state. 

To summarize, both the attenuation and the amplification of light correspond to 
simple analog models that exactly describe the quantum effects of such irreversible 
processes. An absorber is represented by a beam splitter and an amplifier by a 
parametric amplifier. The second mode of the beam splitter represents the absorp- 
tion reservoir, while the additional fluctuation mode of the amplifier describes the 
amplification noise. We have proven |1U21 HU4j the equivalence between these simple 
models and the master equation ()5.1)1 for the Lindblad operators (|5.4J) . i.e. for ther- 
mal reservoirs. We can easily extend jl()21 ll()4j our analog models to phase-sensitive 
Gaussian reservoirs [HE] with the squeezed Lindblad operators 

L 1 =^a + va ] , L 2 = L\, |/i| 2 - \v\ 2 = 1 , (5.18) 

characterized by the complex constants \i and v, see Refs. |1U2[ I104j and Refs. 
cited therein. The transformation (j5.18J) squeezes the thermal Wigner function 
()5.7|) of the quantum noise in one phase-space direction and stretches it in the 
orthogonal direction, indicating that the reservoir is indeed phase sensitive, possibly 
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with reduced fluctuations in one of the quadratures. The fluctuation mode is in a 
state with Gaussian Wigner function and Gaussian density operator Whether 
our simple models can be extended beyond Gaussian reservoirs remains unknown. 



6 Parametric amplifier 

The prime example of an active linear device is the optical parametric amplifier 
168J. Phase-conjugating mirrors |168| and four-wave mixers [168J belong to the 
same category. The quantum optics of parametric amplifiers is studied in Refs. 
^|^E!HIIiaiI32III33llISlin3II2H|, the quantum properties of phase-conjugating 
mirrors are considered in Refs. |2J El 1221 EH CL39J and the quantum effects of four- 
wave mixers are studied in Refs. [8Tl j I157[ I195[ I197j . Active devices require external 
energy sources, often provided by other light beams. These pump beams interact 
with the modes to be amplified in non-linear media [2H H68j , mostly certain crystals 
in the case of parametric amplifiers. The modes are linearly amplified, as long as 
they do not feed back to the pump processes. Here we focus entirely on the regime 
of linear amplification. The quantum physics of pump depletion has been analyzed 



The simplest example for parametric amplification in physics is the playground 
swing. Rocking the legs changes the moment of intertia. Rocking with twice the fun- 
damental frequency of the swing amplifies the oscillation, starting from tiny initial 
movements, a phenomenon called parametric resonance [HE]- The simplest opti- 
cal example of a parametric amplifier is the downconverter |127j . Pump light with 
frequency u p drives two other beams of light, called the signal and the idler, with fre- 
quencies u s + Ui = Up. Assisted by the non-linear medium, some pump photons with 
energy hw p decay into photon pairs with energies fkv s and hu)i. The Hamiltonian 
(J3.23j) describes such a process, the creation of photon pairs. The Hamiltonian also 
accounts for the reverse process where photon pairs become annihilated with their 
energies transferred back to the pump. The relative phase between the signal, idler 
and pump beams decides the direction of the process. Furthermore, momentum con- 
servation requires that the wave vectors k p of the pump light should equal the sum 
of the wave vectors of signal and idler, k s and kj, a condition called phase matching 
168J. Parametric amplifiers with high quantum-noise quality and efficiency tend to 
take pure crystals, good resonators and a number of ingenious experimental tricks. 

6.1 Matrix structure and squeezing 

The Hamiltonian (J3.23|) of the downconverter generates the linear mode transfor- 
mation (|3.1|) with the S matrix ()3.2()j) . The mode operator of the signal, say a±, 
is mixed with the Hermitian conjugate of the idler, a 2 - Simultaneously, the idler 
operator a 2 is mixed with the conjugate of the signal, a\. We would expect the same 
for phase- conjugating mirrors [168 . Let us assume the mode transformation 



in Refs. PBEEI]. 




(6.1) 
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corresponding to the S_ matrix 
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The mode transformation (|6.1|) describes pure amplification, without scattering. 
We require that the mode operators of both the incident and the amplified light are 
proper Bose operators, which results in the quasi- unit arity relation (J3.4j) of the S_ 
matrix. Explicitly, we obtain 
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The general solution of these equations is 

*a/2 ( cosh(9/2) e l ^ + ^/ 2 sinh(6/2) e^-^l 2 
-~ e \ sinl^e^)^^-*)/ 2 cosh(0/2)e- i (* + *)/ 2 

with the real parameters <3>, G, \l/ and A. We express B_ as the product 



B 



JA/2 



; i#/2 






-j*/2 



/ cosh(0/2) sinh(6/2) 
I sinh(6/2) cosh(9/2) 









e -i#/2 



(6.3) 



(6.4) 



(6.5) 



Like in the case of the beam splitter, the net result of the parametric amplifier or of 
the phase-conjugating mirror amounts to four steps. First, the incident beams gain 
a phase $/2, yet in contrast to the beam splitter, this is not a relative phase, but 
an absolute phase shift, because the transformation (|6.1|) acts on a\ and a 2 , not on 
cii and &2. After the phase shift the modes are amplified by the factor cosh(0/2) 
and mixed, with the overlap sinh(0/2), and finally the outgoing modes gain the 
absolute phase \l//2 and the relative phase A. We could include the phases into the 
definitions of the incident and the amplified modes. The hyperbolic mixing (J3.2(Jj) 
would remain as the key feature of the device. We can express the S_ matrix ()3.2()j) 
with ( = B/2 as 



S = R- 
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in terms of the matrix of the 50:50 beam splitter (f4~23jl 
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(6.6) 



V 



l/y/2 l/y/2 

-1/V2 l/y/2 










l/y/2 l/y/2 

-1/V2 l/y/2) 



(6.7) 



In the case when the signal and the idler modes are the linear polarization modes 
of a single light wave the matrix R describes the rotation of the polarization axis by 
7r/4. The parametric amplifier thus processes the rotated modes separately 



VI 



a rl cosh(G/2) + ajj sinh(9/2) 
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a r2 cosh(9/2) - a\. 2 sinh(9/2) . (6.8) 
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We obtain for the quadratures (|2.10j) 



q' rl = qne+®/\ p' rl = p rl e" / 2 , q' r2 = q r2 e^ 2 , p' r2 = p r2 e+ / 2 . (6.9) 

The g ri quadrature is stretched and the p r \ quadrature is squeezed, while preserving 
the Heisenberg commutation relation (|2.11|) . as expected from a canonical trans- 
formation. In the second rotated mode the q r 2 quadrature is squeezed and the p r i 
quadrature is stretched accordingly. Therefore, in the rotated basis, the paramet- 
ric amplifier acts as a perfectly noiseless amplifier or de-amplifier, a squeezer, for 
particular quadrature components. The parametric amplifier may produce squeezed 
light [3l" | 11101 H22j . A reduction by up to — 7dB = 20% in the quadrature variance 
A 2 q compared with the vacuum noise has been observed so far p74"j . 



6.2 Effective Lorentz transformations 

The beam splitter ()4.1)1 generates abstract three-dimensional rotations, expressed 
in terms of the Euler angles $, and \1/ in the decomposition f)4.4j) of the scat- 
tering matrix. To recall a less abstract example, polarizers perform rotations in 
the Poincare sphere, on the quantum Stokes parameters in the Jordan-Schwinger 
representation ()4.6|) . The parametric amplifier ()6.1|) turns out to generate effective 
Lorentz transformations. We define the operators j!98j 

K t = ~ (a{, a 2 ) 1 J = ~ (al^ + a 2 a\ 

K x = ^ (at, a 2 ) a x ( ^{ J = ^ (a x a 2 + a\a\ 

Ky = \ a 2 ) a y J = % - (a x a 2 - a\a\ ^ 

K z = ~ (aj, a 2 ) a z ( jfM = i (ajai - a 2 a|) , (6-10) 

with our notation deviating slightly from Ref. |198 . We see from the Bose com- 
mutation relation that K t = L t + 1/2 and K z = L z — 1/2 where the L operators 
belong to the Jordan-Schwinger representation ()4.6j) . The K operators (|6.1U|) obey 
the commutation relations 

[k x , k y ] = -ik t , [k y , k t ] = ik x , [k t , k x \ = %k y , (e.n) 

while K z commutes with all others. We find 

k 2 t -k 2 x -k 2 y = k z {k z + i), (6.i2) 

the equivalent of the squared angular momentum (|4.8jl . The K operators lie on a 
hyperboloid, see Fig. 01 Since K z commutes with K t , K x and K y , the K opera- 
tors generate transformations with the invariant (J6.12j) . the quantum analog of the 
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Figure 3: Hyperboloid of the parametric amplifier, the counterpart of the Poincare 
sphere of the beam splitter shown in Fig. El The picture schematically illustrates 
the K operators (|6.10p . Depending on the sign of the photon- number difference 
K z + 1/2, the amplifier operates on the upper (positive) or the lower (negative) 
sheet of the hyperboloid. Amplification processes draw trajectories on one of the 
sheets, consisting of the Lorentz transformations ()6.14|) . 



squared space-time distance [HZ] in 2+1 dimensions. In other words, the K operators 
generate effective Lorentz transformations. In fact, we find 
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two Lorentz transformations jHZ] with effective velocities v x /c = tanhO and v y /c = 
tanhfi and one rotation with angle $. We obtain for the mode operators 
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(6.14) 



These are all transformations that belong to the class (|6.1|) . Therefore, parametric 
amplifiers generate effective Lorentz transformations. The measured quadrature- 
noise reduction (Mj of 7dB (20%) gives, according to Eq. (|6.9|) . a 9 parameter of 
— log 0.2, which corresponds to an effective velocity of v/c = tanh G = 0.92. Finally, 
we employ the Lorentz generators (j6.10|) to represent the evolution operator of the 
parametric amplifier with the matrix (j6.5|) . We use our results f|6.14|) and get 



B = exp(-i&K t ) exp(-iQKy) exp(-i^K t ) exp(-iAK z 



(6.15) 



Passive optical instruments such as beam splitters preserve the total number of 
photons ()3.7|) . they conserve energy. The parametric amplifier is subject to a con- 
servation law as well. The generators of the evolution operator (J6.15j) commute with 
K z . Consequently, the photon- number difference is conserved, 



i't; 



i't; 



a\(ii 



(6.16) 



Photons are emitted in pairs, one in the signal and one in the idler beam, as we 
would expect from the Hamiltonian ()3.23j) . 



6.3 Quantum correlations 

Conservation laws manifest themselves in correlations. For example, picture a pair 
of two particles with spin, say two polarized photons. The particles are produced 
at one spot and then they move away from each other. Suppose that the total 
spin of the two particles has been zero initially and that spin is conserved. Now, 
if we measure the spin of the first particle, the second one must have the opposite 
spin, regardless which type of polarization we are probing, linear or circular, and 
regardless how far apart the particles are. Such long-ranging correlations exist due 
to conservation laws, they do not violate the relativistic causality, because they do 
not cause each other, but rather have a cause in common, and they are completely 
classical. Quantum mechanics adds a subtlety, a decisive yet quantitatively subtle 
feature. For quantum particles the outcomes of the spin measurements may be 
probabilistic, depending on their state, but the spins of the partners are correlated. 
For example, consider the singlet state of two photons in two light beams, with two 
orthogonal polarization modes (+) and (— ) each, 

W = -4(l+>i|-> a -(|->i|+> 2 ). ( 6 - 17 ) 
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or, written as Fock states, 

l^) = ^ f (|l)i + |l)2--(|l)i-|l) 2+ ). (6.18) 

The polarization of photon 1 is completely random, (+) or (— ) with 50% chance, 
but the polarization of the second photon is always correlated to the first one. 
The quantum state (|6.17|) is entangled — the vector \ip) of the total state does 
not factorize into the state vectors of the subsystems. Now, suppose we measure 
the polarization of the first photon, with respect to a given axis, and we turn the 
polarizer of the second by an angle. In this case, the measurement results are not 
perfectly correlated anymore. However, we can quantify the correlation degree by 
calculating the difference D between the statistical frequency of coincidences and the 
frequency of detecting the opposite spin, after many repeated runs of the experiment 
for each setting of the polarizer angle. Any spin measurement is characterized by 
the scalar product of a three-dimensional unity vector a with the vector of the Pauli 
matrices a = (a x ,a y ,a z ) T of Eq. (|4.5j) . where a describes the polarizer setting on 
the Poincare sphere. We obtain for the singlet state (jbM7j) 

D(a u a 2 ) = (ip | (a : • a x ) (a 2 • <j 2 ) I V>) = ~ a i ' a 2 • (6-19) 

It turns out fDJ 120] that D violates certain inequalities for local hidden variables, 
based on classical statistics, called Bell's inequalities after their discoverer [Tj5]. Sup- 
pose that the particles are classical by nature, but are engaged in a quantum con- 
spiracy. Before they are separated they agree on the outcome of spin measurements 
denoted in some hypothetical hidden variables that determine the measurement re- 
sults. Although the observer is unaware of the plot, he can put constraints on the 
correlator, for example 02] 

D( ai , a 2 ) + D( ai , a' 2 ) + a 2 ) - D(a£, a' 2 ) < 2 , (6.20) 

assuming that the particles cannot communicate with each other after they have 
been separated. The spin correlations (jSHSj) of the quantum state (|6.17|) violate the 
Bell inequality ()6.20|) by maximally a factor of \/2, as one may verify Quantum par- 
ticles can be slightly stronger correlated than classical statistics allows. Violations 
of Bell's inequalities have been convincingly demonstrated in several experiments 
since the major breakthrough in efficiency due to Aspect 0, see e.g. Refs. |174|ll90j . 
The parametric downconverter |127| Q68 a is the central tool in most modern tests 
of genuine quantum correlations. Apart from probing the foundations of quantum 
mechanics, quantum correlations are beginning to play a decisive role in tasks be- 
yond the capabilities of classical physics such as quantum cryptography [SHI EHj 
and in other forms of quantum communication [3U1 I135j . Parametric amplifiers are 
frequently applied in this wider field of entanglement engineering. 

The debate on quantum correlations, on the "spooky action at a distance" , that 
preceded any applications of entanglement, dates back to the 1935 paper [5Tj "Can 
Quantum- Mechanical Description of Physical Reality Be Considered Complete" by 
Einstein, Podolsky and Rosen. Schrodinger |163j introduced the term entanglement 
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to physics, Verschrdnkung, a German term used by cabinet makers for dovetailing, 
quite a fitting term for the quantum correlations in the natural world. The paper 
|o"T] by Einstein, Podolsky and Rosen has spawned a large literature |129| and has 
inspired a series of intriguing experiments, see for example Refs. [HI I2H I2H HUB H4Ut 
I141[ 11741 119U| . Interestingly, the parametric amplifier represents not only the device 
of choice for many modern applications of quantum correlations, but the amplifier 
is also able to reproduce the original argument by Einstein, Podolsky and Rosen for 
the quadratures q and p |158j in the experiment |14()| I141j . 

Suppose the signal and idler modes are in the vacuum state initially. An opti- 
cal parametric amplifier is just an amplifier after all. So, if the initial amplitude 
is zero the amplified amplitude remains zero. However, amplitude fluctuations get 
amplified. Remember the classic example of parametric resonance [HE], the child on 
the playground swing where rocking at twice the frequency of the swing amplifies 
tiny initial movements. In the case of the optical parametric amplifier, the vac- 
uum fluctuations are amplified. The resulting quantum state, called the two-mode 
squeezed vacuum, comes close to the original state considered by Einstein, Podolsky 
and Rosen in their 1935 paper jHI]. We calculate its Wigner function. Initially, the 
two modes are in the vacuum state with the Wigner function |ll()j 



W(q 1 ,p 1 ,q 2 ,p 2 ) = — exp {-q\ - p\ - q 2 - p 2 2 ) 
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(6.21) 



We recall that simple optical instruments transform the Wigner function as if W were 
a classical probability distribution of the quadratures. For simplicity, we assume that 
the S_ matrix ()6.2|) of the amplifier is the real matrix ()3.20|) . We obtain from the 
transformation rule ()3.30|) the Wigner function of the two-mode squeezed vacuum 
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Consequently, 



W — — exp 

TT 2 



(gl ~ g2) 

2e- 2 f 



(gi + <fe) 2 _ (gi +P2) 2 _ (P1-P2) 2 ' 
2e 2 < 2e- 2 ^ 2e 2 < 



(6.24) 



The Wigner function ()6.24|) indicates strong quadrature correlations when \(\ is 
large. Suppose we measure the q quadrature of the first beam using homodyne de- 
tection |llUj and we find the value q±. In this case a quadrature measurement of the 
second beam, some distance away from the first one, gives with high probability q\ 
as well. The first measurement has prepared the second beam in a state that closely 
resembles a q quadrature eigenstate. Now, if we measure the p quadratures, with 
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the result p\ for the first beam, the second mode produces —p\ with high probabil- 
ity. The q and p eigenstates cannot coexist, because of the canonical commutation 
relations (J2.11j) . Therefore, depending on the measurements on the first beam, the 
second one is prepared in states that are mutually incompatible and that cannot 
coexist as "elements of reality" [HI]. Of course, Einstein, Podolsky and Rosen did 
not mention quantum-optical quadratures in their paper [SI], but rather position 
and momentum in real space, yet this does not fundamentally alter the argument. 
However, a hypothetical quantum conspiracy could explain the behavior Einstein, 
Podolsky and Rosen found so puzzling. The light beams could have originally agreed 
on with which probability they are going to produce all linear combinations of q and 
p quadratures. The Wigner function ()6.24|) reveals the plot. The function plays the 
role of the non-negative probability distribution for some hidden variables trigger- 
ing all quadrature detection results. In our case, the quadratures themselves are 
the hidden variables qi, q2, pi and p 2 . Ironically, the paradox that led to quantita- 
tive measures of quantum nonlocality can be explained classically. It would take a 
negative Wigner function |lU7j to violate Bell's inequalities in quadrature space, as 
Bell has pointed out [21], although with flawed mathematics in his counterexample 
|82j . However, the two-mode squeezed vacuum state may well violate Bell's inequal- 
ities for other observables than the quadratures [H], because the quasiprobability 
distributions in other variables [10SJ I109| I182j may be negative. 

The parametric amplifier generates strong quantum correlations between the 
emerging light beams. In fact, one can argue [T3l IT^j that the correlations are the 
strongest possible for two equal oscillator modes with given average energy. To 
understand why, we calculate the reduced Wigner function of any of the modes, say 
mode 1. We obtain after Gaussian integration 

/+oo P+OO 
/ W'(q 1 ,p 1 ,q 2 ,p 2 )dq 2 dp2 = W th (q u p 1 ) , (6.25) 
-oo J — oo 

where W t h denotes the Wigner function ()5.7|) of the thermal state ()5.8|) with the 
temperature 

T= —-^— . (6.26) 
2/c B lncothC 

The reduced Wigner function describes the reduced quantum state of one of the 
correlated modes, averaged over the second mode. As we know from statistical 
physics [99J, the thermal state is the state of maximal entropy with given energy, 
in the canonical Gibbs ensemble, the state of maximal disorder and minimal infor- 
mation. On the other hand, the total two-mode squeezed vacuum is a pure state. 
Consequently, this state stretches mostly across the correlations between the two 
entangled modes. When reduced to one mode, the subsystem contains as little in- 
formation as possible. Therefore, the two-mode squeezed vacuum is indeed the most 
strongly correlated quantum states of two equal harmonic oscillators with given en- 
ergy |T2| E] • The higher this energy the higher is the effective temperature (J6.26|) 
and the stronger are the quantum correlations. 

Since the parametric amplifier conserves the photon-number difference f!6.16|) 
between signal and idler, the quantum correlations of the two-mode squeezed vac- 
uum are particularly transparent in terms of photon- number states, i.e. in the Fock 
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representation. We utilize the disentangling theorem I147j 



exp(2i(K y ) = exp(a\al tanh£) exp(— 2K t In cosh £) exp(— a\a,2 tanh£) , (6.27) 



and get for the two-mode squeezed vacuum state 



\if>) = exp(2«CA v ) |0)i|0) 2 = exp(ala|tanhC) exp(-2K 4 IncoshC) |0> 1 |0> 2 



We see that consists entirely of photon pairs |n)i|n)2. We also see that the 
reduced density operator tr2{|'0}('0|} is indeed the thermal state (jf>.8j) with tem- 
perature (|6.26|) . In spontaneous parametric downconversion |127j the squeezing 
parameter ( is small. In this regime, ( 2 gives the probability for producing a single 
pair of photons, as we see from a first-order perturbation theory with the Hamilto- 
nian (J3.23|) . and ( 2n gives the probability for n independently produced pairs. So, 
according to the result (|6.28|) . the two- mode squeezed vacuum seems to consist of 
indistinguishable photon pairs, in the regime of weak amplification. Spontaneous 
downconversion |127j has been frequently used for generating single pairs of pho- 
tons, because the probability for multiple pairs is very low, ( 2n compared with ( 2 . 
High-efficiency single-photodetectors operate like Geiger counters, where a single 
photon triggers an electron- hole avalanche |llUj . Since the initial photons are just 
triggers, the detectors cannot discriminate between one or more detected photons in 
each detection click. Therefore, it is an experimental advantage to guarantee that 
mostly single pairs occur, which of course limits the rate of pair production. 

The interference experiments with single photons mentioned in Sec. 4.3 have been 
performed with photon pairs generated in spontaneous parametric downconversion 
[127J. Here the quantum state (|6.28j) of light is essentially 



In such experiments only those experimental runs count where photons are counted, 
the time when the detectors are not firing is ignored, which reduces the quantum 
state to the photon pair |1) 1 |1) 2 . Postselection disentangles the two-mode squeezed 
vacuum. We argued in Sec. 4.3 that the interference of the photon pair |1) 1 |1) 2 at 
a 50:50 beam splitter generates the entangled state (|4.24p . Without postselection, 
however, this state is the disentangled product of two single-mode squeezed vacua, 
as we see from the factorization (|6.6|) of the S_ matrix. The notion of entanglement 
is to some extent relative. 

Finally, we explain briefly how the singlet state (j6.18j) is generated in spontaneous 
parametric downconversion |127j . The spin singlet ()6.18|) refers to the polarization 
state of two light beams with two polarization modes each, four modes in total with 
the mode operators Oi + , ai_, a 2+ and a 2 _. The singlet is completely unpolarized, 
but contains polarization correlations, Bell correlations i n fact. Therefore, the 





(6.28) 



|V>)«|o>i|o) 2 + c|i>i|i>2- 



(6.29) 



36 



creation process of the singlet must be polarization invariant. The Hamiltonian of 
parametric downconversion is quadratic in the mode operators of signal and idler, as 
long as the pump can be regarded as a classical parameter. Consequently, we need a 
quadratic form of the annihilation operators a 1+ , di_, a 2 + and a^- that is invariant 
under polarization transformations |lUlj . The rotations on the Poincare sphere ()4.9|) 
correspond to the mode transformations (|4.10J) of the beam-splitter type (|4.1J) of the 
polarization modes (ai+,ai_) T and (0-2+, a^-) 7 ', with identical SU(2) matrices j4"K] . 
Written in matrix form 

a' = Ba, a=( ° 1+ ° 2+ V B' 1 = B\ detB = 1 . (6.30) 
- - V. a i- fl 2- J ~ ~ 

The determinant of the a matrix has the desired property, because deta' = det-B deta. 
Consequently, we arrive at the effective Hamiltonian 

H = ie iip ( deta f - ie'^C deta (6.31) 

with real squeezing parameter ( and real phase (p. We assume that both signal and 
idler have initially been in the vacuum state. We use the result ()6.28|) and get 

|V) = expH^)|0) 1+ |0) 1 _|0) 2+ |0) 2+ = |^) + + |V)_ ) 
1 00 

\iph = ^h7E( ±e_inanh n |n)i±|n)2T, (6.32) 

* n=0 

which, after postselecting for single photon pairs, gives the Bell state (j6.18|) in the 
regime of spontaneous parametric downconversion. We describe briefly how the 
downconverter can realize the Hamiltonian (|6.31j) . 

Downconversion is a nonlinear optical process fZG\ I168j where a pump beam 
E p of frequency u p generates two waves, the signal E s and the idler E iy of fre- 
quencies u s and Ui with uo p = uo s + Wj. This three-wave mixing [16*5| requires an 
appropriately anisotropic medium, a crystal, because the interaction Hamiltonian 
of the three fields y^ 1 E v E s Ei is not invariant under parity transformations where 
the fields change sign. In an anisotropic medium light propagates as ordinary or 
as extraordinary waves with opposite polarizations |27j . There are two principal 
types of downconverters. In type I [IBS] the modes E s and are both ordinary 
or extraordinary waves. In type II downconversion |168j either E s or E{ is ordinary 
and the other wave is extraordinary, and hence the two have opposite polariza- 
tions. The downconverted light leaves the crystal in two cones that display the 
conservation laws involved — the conservation of energy, uo p = uo s + uii, and the con- 
servation of the wave vectors (phase matching condition |168j ). see Fig. 0] Where 
the cones intersect the polarization state is undecided. According to our theory, this 
polarization- invariant state of light is the state ()6.32|) . The two intersection lines of 
the emerging light cones carry polarization Bell states [01], see Fig.EJ Alternatively, 
one can employ two subsequent type I crystals with orthogonal optical axes |H21- 
The crystals downconvert pump light into entangled beams in all directions. All 
this shows that the parametric amplifier is a superb device for generating quantum 
correlations of light. 
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Figure 4: Schematic of the downconversion process (with type-II phase matching). 
An ultraviolet photon incident on a nonlinear crystal can sometimes spontaneously 
split into two correlated infrared photons. These photons are emitted on opposite 
sides of the pump beam, along two cones, one of which has horizontal polarization, 
the other of which has vertical polarization. Photon pairs emitted along the inter- 
sections of the cones are entangled in polarization — each photon is individually 
unpolarized, and yet the photons necessarily have perpendicular polarizations, no 
matter how far apart they are. (Courtesy of Paul Kwiat.) 




Figure 5: Entangled light. Colorized infrared photograph of the measured downcon- 
version output jUTj schematically shown in Fig. |3J The green rings correspond to 
photons with roughly equal energies (half that of the parent pump photon). Where 
the rings intersect the polarization state of the light is entangled. (Courtesy of 
Anton Zeilinger, Institut fiir Experimentalphysik, Universitat Wien.) 
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Figure 6: Black hole. Gravity modifies the metric of space and time. In flat space, 
far away from the gravitating object, the light cones stand upright in a space-time 
diagram. Gravity is tilting the light cones, bending light. At the event horizon 
all future light cones point inwards. Light cannot escape anymore, and so cannot 
anything else. (Picture by Enrique Arilla.) 

7 Optical black hole 

Black holes are optical instruments as well, although they tend to operate on grander 
scales than ordinary laboratory equipment. First of all, black holes are gravita- 
tional lenses |162j - the gravitational fields surrounding them focus light from 
other sources and may cause a multitude of optical illusions. Black holes share this 
feature with other gravitating bodies such as entire galaxies, stars or even individual 
planets. The distinctive feature of a black hole is the event horizon |131j . The hori- 
zon, surrounding the hole, is a place of no return. Beyond the horizon, everything 
falls into the hole, no matter how fast it moves, even light. Here gravity has tilted 
the future light cones such that they all point inwards, see Fig. |U1 

As we discuss in this section, the horizon turns the black hole into an active 
optical instrument. The hole does not merely act as a gravitational lens, it emits 
radiation as well, Hawking radiation [23 HUH E3 EE] • The hole appears as a thermal 
black body with a temperature that depends on the gravity at the horizon (and hence 
|131j on the horizon's area and the hole's mass). The characteristic wavelength at 
the peak of the Planck spectrum is comparable to the radius of the horizon [23 EH] , 
more than 1km for the known solar-mass or larger black holes, which makes Hawking 
radiation far too cold to be observable against the cosmic microwave background of 
about 1mm wavelength. 
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On the other hand, since black holes are optical instruments after all, why not 
use optical equipment to make one? Such table-top black holes would be much 
smaller than the astronomical ones. Since the scale of the hole determines the 
temperature of the Hawking radiation, one would expect them to radiate stronger 
than their larger archetypes in space. Most of the proposed artificial black holes [136J 
are based on a simple idea |181| 1184) : Consider a moving medium with spatially 
varying flow speed, say water going down the drain of a bathtub. Waves in the 
fluid, sound waves |181| I184j . surface waves [165] or light waves |3H H54] are 
trapped when the flow exceeds the speed of the wave in the medium, the speed of 
sound or the effective speed of light, for example. To turn this simple analogy into 
measurements of emergent Hawking radiation will take extraordinary media though 
- Bose-Einstein condensates of alkali vapours [HIJ 1154) for sonic holes 1117] , 
two phases of superfluid Helium-3 |18fi| for ripple holes |165j and slow-light media 
|116t 1125] for optical holes |112] 1114"! I115j . Note that the possibly realistic version 
of the optical hole |114t I115j differs from ordinary black holes. It does not involve 
a moving medium and belongs to a different class of quantum catastrophes |115j 
than the Hawking effect [2EI ESI EH EEJ- Much in this field is still in flux and 
is connected to other areas of physics outside of the scope of this article. Here 
we discuss an optical Gedankenexperiment, not a realistic proposal, to show how 
Hawking radiation emerges in principle. 



7.1 Light in moving media 

Consider a dielectric medium with an incredibly strong electric permittivity e ^> 1 
and with unity magnetic permeability, for simplicity. Such a medium would signifi- 
cantly reduce the phase velocity of light, c', because 

r 2 

c ,2 = -. (7.1) 

e 

Note that the slow light |ll(j| 1125] demonstrated in recent experiments [UHl 1121] I153j 
represents a different case where the group velocity of light is reduced and not the 
phase velocity. Consider for simplicity a one-dimensional medium moving with the 
spatially varying flow speed u(x). The Lagrangian of the electromagnetic field is the 
Lorentz scalar JHH II 13j 

jgf = I (ED - BH) . (7.2) 

At any place in the moving medium we can imagine a co- moving frame of coordinates 
where the medium is at rest, denoted by primes. The Lagrangian is invariant, 

Jgf = - (E'D' - B'H') . (7.3) 

In a local co-moving frame, the D' and H' fields are connected to the E' and B' 
fields via the constitutive equations [HJ HUUj 

D' = £ eE', H' = e c 2 B' . (7.4) 
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We represent the E' and B' fields in terms (12. 2 j) of the vector potential A, a scalar 
in 1+1 dimensions. For non-relativistic velocities of both the medium and the light 
we obtain from the local Galilei transformations to the laboratory frame 

& t = dt + ud x , d' x = d x . (7.5) 

In this way we arrive at the Lagrangian in the laboratory frame 

J8f = £ -f [{d t Af + 2u(d t A)(d x A) + {u 2 - d 2 ){d x A) 2 ) (7.6) 

that generates the Euler-Lagrange equation [§7] 

d t e(d t + ud x )A + d x e[ud t + (u 2 - c' 2 )}d x A = . (7.7) 

This wave equation describes the propagation of classical light in moving media in 
one spatial dimension. 

7.2 Space-time geometry 

Remarkably, light turns out to experience the moving medium as an emergent space- 
time geometry [H21 DHH Hl^ I137| I151| I152j . To understand why, consider the sim- 
plest case, our one-dimensional medium with large e. We introduce the relativistic 
notation 

x° = t, x 1 = x, d = d t , di = d x (7.8) 

and adopt Einstein's summation convention over repeated indices. In this way we 
write the Lagrangian ()7.6|) as 



J? = jr(d,A)(d u A) (7.9) 



in terms of the matrix 

which gives the Euler-Lagrange equation 

a^ra v A = Q. (7.ii) 

Suppose that e is constant while u may vary in space. We introduce the metric 
with the determinant [HZ] 

9 = -c' 2 (7.13) 



and the inverse 

1_ f 1 u 
c' 2 V u -c' 2 + u 



r = -A : ) ■ (7.i4) 
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As long as e remains constant we can write the wave equation (j7.11|) as [97 



The Dp denote the covariant derivatives and D^D^ gives the D'Alembert operator 
in general coordinates with the metric ([7.12)1 . Equation (|7.15|) describes a mass- less 
field, light, in general relativity. Therefore, the electromagnetic field perceives the 
moving medium as the space-time geometry ([7.12)) . Light-rays follow zero-geodesic 
lines measured with respect to the metric ()7.12[) |lllj . In the history of optics, the 
oldest known paradigm of geometric guidance is Fermat's Principle |27j : Light rays 
follow the shortest optical paths where the path length is measured with respect to 
the refractive index. 

A dielectric medium acts on light, in focusing or scattering light, but, according 
to the principle of actio et reactio, light also acts on the medium, in exerting optical 
forces. We have ignored these forces here, assuming that the medium is very heavy. 
The light forces turn out to follow geometric ideas, too. The medium perceives the 
electromagnetic field as a space-time geometry as well. Light and dielectric matter 
see each other as geometries |113j . 

Finally, the idea that media modify the geometry of space and time can be turned 
around, at least speculatively. Can one interpret space-time curvature, i.e. gravity, 
as generated by a hypothetical fluid, as an emergent phenomenon in condensed- 
matter physics |186j ? Sir Isaac Newton speculated in 1675 that gravity acts 
through "an aetherial medium, much of the same constitution with air, only the 
vibrations far more swift and minute" , one of Sir Isaac's departures to hypotheses. 
Much later, in 1968, Andrei Sakharov, father of the Soviet hydrogen bomb, dissident 
and Nobel Peace Price laureate, noticed |16Uj that Einstein's general relativity could 
appear as the effect of fluctuations of the quantum vacuum, the modern form of the 
aether. 

7.3 Horizon 

The effective metric of light in moving media |6^ I113j does not cover all possi- 
ble space-time manifolds, but the generated geometries are rich enough to include 
horizons. The horizon is the place where the flow speed u exceeds the speed of 
light in the medium, c'. Here, according to the metric ()7.12)) . the measure of time, 
c' 2 — u 2 , is zero — time appears to stand still. Light propagating against the cur- 
rent freezes. In fact, black holes have been termed "frozen stars" jHS]. Immediately 
before the horizon, any outgoing light struggles to escape. The closer the distance 
to the horizon, the more optical cycles it takes to move forward. In turn, the wave 
length of light shrinks to zero at the horizon. On the other side of the horizon, no 
counter-propagating light is able to escape anymore. 

Let us turn these words into formulae. Assume for simplicity that e is constant. 
In this case the general solution of the wave equation ()7.11)) is 




(7.15) 





(7.16) 
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Note that also the complex conjugate, A*, is a solution, because the wave equation 
(|7.11|) is real. Suppose, without loss of generality, that the medium flows from the 
right to the left such that u is negative. Hence A + (r + ) describes counter-propagating 
and A_(r_) co-propagating wavepackets. The t± are the null coordinates in the 
global medium frame. (Note that in more than one spatial dimension such a frame 
does not exist in general.) Assume that the velocity u exceeds the effective speed of 
light at one point, say x — 0, such that 

u ~ — c' + ax near x = . (7-17) 

The positive constant a describes the velocity gradient at the horizon, the analogue 
of the surface gravity at the event horizon of a gravitational hole. In the vicinity of 
the horizon we get 

r + ~t- Hx/X±oo) (7.18) 
a 

with the integration constants o;±oo- Therefore, to escape from the horizon takes an 
exponentially long time. Assuming monochromatic light, where A + = exp(— iur + ), 
we see that the phase <p of the light wave is logarithmic in x/x^. Consequently, the 
wavelength 27r/(d x ip) is indeed proportional to the distance from the horizon, 

A-Aoo— . (7.19) 

We also see that the r + coordinate distinguishes the left and the right side of the 
horizon. The constants x± OQ must be positive on the right and negative on the left 
side, for not running into problems with complex time. The horizon cuts space into 
two disconnected parts. 



7.4 Hawking radiation 

Consider the quantum theory of light for our artificial black hole, the moving 
medium. As usual [23 EHj- the quantum field A consists of a set of modes with 
the mode functions Ak and the Bose annihilation operators a&, see Eq. ()2.f)jl . In the 
moving medium, the mode functions are subject to the wave equation (|7.11j) . The 
modes are normalized according to a scalar product that ought to stay invariant in 
time and that should approach the scalar product (|2.4J) in the limiting case when 
the medium is at rest. We put 

(A 1 ,A 2 ) = (A\f v d u A 2 -A 2 f v d u A^dx 

= Y J (A\{d t + ud x )A 2 - A 2 {d t + ud x )A\)edx. (7.20) 

The scalar product ()7.20|) is the desired one, because it agrees with Eq. (12. 4j) when 
u vanishes and the product is indeed a constant in time, because of the conservation 
law 

<9 M (Al p v d v A 2 - A 2 rduA\) = . (7.21) 
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After these general remarks on quantum light in moving media we return to the 
optical black hole. 

Consider the monochromatic modes with frequency u that are propagating against 
the current, the mode function Ar on the right side of the horizon and the function 
A L on the left side, 

Ar = A R Q(x) exp(-iujT + ) , A L = A L Q(-x) exp(zcjr + ) . (7.22) 

The unity step functions 9(±x) indicate that A R is only supported on the right side, 
whereas Al exists solely on the left side. Note that the mode Al on the left side 
beyond the horizon oscillates with the negative frequency — u, in order to represent 
a proper positive-norm mode (|2.7jl . In fact, we see from 

i(d t + ud x )A L = -uo -4-: A L (7.23) 
u + c 

that the normalization integral ()7.20|) is positive for negative u + c', as is the case 
beyond the horizon. Negative frequencies u> correspond to negative energies hu. 
Consequently, the artificial black hole does not have a natural ground state, at least 
in our simplified model where the permittivity e and the flow u are given and fixed, 
ignoring the backaction of the light onto the medium. The vacuum state, the state of 
zero photons, traditionally the ground state, may thus depend on the history of the 
horizon, on the catastrophic event creating the black hole. In space, the gravitational 
collapse creates the black hole. In our laboratory analogue, the medium must have 
been accelerated in the past to form a stationary flow with a horizon. Consider the 
quantum theory of light in the Heisenberg picture where the field operators evolve 
and the quantum state of light remains constant. Trace the mode functions back 
into their pre-horizon history, in order to find the decisive signature of the modes 
that distinguishes the initial quantum vacuum, and hence the final vacuum, too. 
Complex analysis provides us with elegant theoretical tools for this purpose ISO] • 
Before the horizon has been formed, the vacuum modes stretched over the entire 
medium. Counter-propagating waves consisted of superpositions of plane waves with 
positive wavenumbers, 

r+oo ^ 

A (x, t) ~ / A (k) exp(ikx - iut) dk , t -> -oo . (7.24) 
Jo 

Consider the analytic continuation of the vacuum modes to complex x values. The 
A are analytic on the upper half plane, because here exp(ikx) decays for positive 
wavenumbers k, such that A cannot develop singularities. Now, imagine that the 
original modes evolve during the formation of the horizon where the medium is 
gently accelerated to exceed the effective speed of light beyond one point. A smooth 
process cannot fundamentally alter the analytic properties of waves. Consequently, 
we can employ the analyticity of the mode functions as a marker to distinguish the 
vacuum state of light. The modes Ar and Al are clearly non-analytic, because they 
lack support on either one of the horizon sides. Yet we can glue Ar and A* L together 
to form analytic wave functions Close to the horizon, we have 

exp(-^r + ) ~ (x/x +00 ) iuj/a exp(-iut) . (7.25) 
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Circumventing the horizon on the upper half plane, we use the property (— z) lfl = 
e _7r/i z* M , and find the combinations 



A x = A R cosh C + A* L sinh C, A 2 = A* R sinh ( + A L cosh ( (7.26) 

that are analytic on the upper half plane when 

tanhC = e-™ /a . (7.27) 

The vacuum modes Ai and A 2 are normalized according to Eq. (|2.7|) with the scalar 
product (J7.20J) . taking into account that both Ar and A^ are normalized and orthog- 
onal, and that the norm of A* is the negative norm of A. The mode transformations 
(fT2fi|) imply that 



or \ / cosli^ — sinh£ \ f o-i 
a[ ) \ — sinh ( cosh ( 



(7.28) 



The black hole acts as a parametric amplifier, transforming the incident vacuum into 
the two- mode squeezed vacuum state (jd28j) . the entangled Einstein- Podolsky- Rosen 
state [HI] • The steeper the velocity gradient a the stronger is the entanglement and 
hence the energy of the two-mode squeezed vacuum. The horizon generates photon 
pairs, both propagating against the current, but one escaping on the right side of 
the horizon and the other drifting to the left. Somehow the flow must provide 
the energy for the photon production that, in turn, will cause friction. In space 
the gravitational field surrounding the black hole, created by the mass of the 
collapsed object, should provide the energy for the emerging Hawking radiation. 
The black hole is evaporating and is getting smaller. All black holes are similar, 
apart from their scales |131j . Therefore, smaller holes generate larger gradients of 
the gravitational potential. Consequently |67j . the evaporation is accelerating, until 
the black hole explodes in a grand fiery finale. 

If we confine ourselves to the right side of the horizon, our natural place in the 
case of gravitational black holes, we perceive the two-mode squeezed vacuum f!6.28|) 
as the thermal state ()5.8|) with the temperature f)6.26|) . the Hawking temperature 

T = ^- . (7.29) 

The temperature depends only on the velocity gradient a, the analogue of the surface 
gravity of the hole. Consequently, the Hawking temperature (j7.29|) is the same for 
the entire spectrum, which is highly non-trivial. The black hole appears as a black- 
body radiator in thermal equilibrium [T71 [TBI ESI EH] , a mysterious result [BTj . 



7.5 Grey-body factor 

We have shown that the black hole acts as an active optical instruments, as a para- 
metric amplifier. So far, we have ignored the passive optical properties of the hole, 
the analogue of the gravitational lens |162j . Here we sketch how these properties 
are taken into account, i.e., how to include the light scattering of the optical black 
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hole. In our simple one- dimensional example, scattering occurs only when the per- 
mittivity e of the medium varies in space, because the solution (J7.16)) for constant 
e maintains co- and counter-propagating light as the separate waves A + (t + ) and 
A_(t_). When e varies, we represent the monochromatic modes on the two sides of 
the horizon as 

Ar = A +R exp(-iur + ) + A_ k exp(-i^r_) , 

A L = A +L exp(+iur + ) + A_ k exp(— iujtJ) , (7.30) 

where the A + r, A + l and A_ k are slowly varying functions of x. Similar to Sec. 2.3, 
we construct the transfer matrix (j2.16|) in order to connect A + r and A_j, between 
the various points on the right side of the horizon. Also in the case of moving media, 
the transfer matrix exhibits the structure (j2.17j) with the determinant (j2.18|) . The 
latter property is a consequence of the conservation law (|7.21j) that implies for 
monochromatic modes 

d x (t{eu 2 - c 2 ){A* d x A - Ad x A*) + 2ueu\A\ 2 ^) =0. (7.31) 

Substituting the structure (|7.30j) with the transfer-matrix connection (j2.16J) gives 
the differential equation for the determinant ()2.18|) . 

Similar to Sec. 2.4, we use the transfer matrix to define the in and out modes, 
but starting here with the outgoing modes, 

A° R ut = u R (x) e~ iuJt , A°^ = u D {x) e~ iuJt , (7.32) 

where A ^ describes an outgoing counter-propagating wave, whereas refers to 
a co-propagating wave. The A° R nt wave originates at the horizon, propagating to the 
right, whereas A ^ propagates to the left at +oo and crosses the horizon without 
anything unusual happening. We assume that at +oo the effective speed of light, d 
and the flow u approach constants such that light approaches plane waves with the 
wavenumbers 

kn — — -. r -. r , kjj = — r -. r . (7.33) 

c(+oo) + u[+oo) " c'(+oo) — m(+oo) 
In this case ur obeys the asymptotics 

u nW~ A *{ a * e ^llZ-^ : IZ^L (7 ' 34) 

and ud satisfies the corresponding asymptotics similar to the complex conjugate 
of Eq. (j2.21j) . On the other side of the horizon, we define an outgoing counter- 
propagating wave by 

A L = u L (x) , u L (x) ~ A L (x/x^) 1 ^ : x -> -0 . (7.35) 

We find that Af^ is orthogonal to A R ut with respect to the scalar product (|7.20j) 
on the right side of the horizon. Furthermore, Aff 1 is automatically orthogonal to 
^out on l e ft side, because here the two waves oscillate at opposite frequencies. 
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Consequently, the A^, A^, establish a set suitable for the mode expansion 
()2.6|) of quantum light, the set of the outgoing modes. The ingoing modes are the two 
analytic vacuum modes A x and A 2 of Eq. (J7.26|) . and a mode A 1 ^ that is connected 
to A^ and A™ 1 by the scattering matrix B_ of Eq. ()2.26|) on the right side of the 
horizon. We thus arrive at the transformations 



Consequently, the performance of the black hole as a passive optical instrument is 
entirely due to the light scattering on the right side of the horizon. The scattering 
matrix describes the fraction ()2.26|) of the generated Hawking radiation that reaches 
the detector at x — +00. The rest of the generated light is scattered into co- 
propagating waves across the horizon and lost. The black-body radiation of the 
hole is thus moderated by the grey-body factor ct _1 \/|a| 2 — \b\ 2 . 

8 Summary 

We have gone a long way from studying a piece of glass to the quantum physics of 
black holes. Objects as diverse as beam splitters, multiports, interferometers, fibre 
couplers, polarizers, gravitational lenses, parametric amplifiers, phase-conjugating 
mirrors and black holes have something in common — they act on light as sim- 
ple optical instruments, i.e. as linear optical networks. A linear optical network 
turns a set of incident light modes into an equal number of outgoing modes by a 
linear transformation. If the transformation does not mix annihilation and creation 
operators the instrument is passive and it does conserve the total number of pho- 
tons. Otherwise, the instrument is active and requires an energy source (or drain). 
Such optical networks are reversible devices where, in principle, all the outgoing 
modes can be reversed to restore the incident modes. However, we can also model 
irreversible devices such as absorbers or amplifiers by fictitious beam splitters or 
parametric amplifiers, respectively, where we keep track of only one of the outgoing 
modes. The quantum noise of the unobserved modes accounts for the absorption or 
amplification noise and the mode itself represents the absorption or amplification 
reservoir. 

Our starting point has been the quantum optics of a dielectric slab, say a piece 
of glass. This simple example indicates how to develop a quantum theory of light in 
dielectric media and it shows how the mode transformations of passive optical instru- 
ments come about. We have developed the theoretical tools to analyze the quantum 
physics of linear optical networks and we have applied them to a few characteris- 
tic situations, to the splitting and interference of photons and to the manifestation 
of quantum correlations in parametric downconversion. We have sketched how to 
describe irreversible processes in quantum mechanics and how to apply this theory 
to determine effective models for absorbers and amplifiers. Finally, we returned to 
the starting point, revisiting the quantum optics in dielectrics. We showed that the 
creation of a horizon turns a passive medium into an active device, similar to the 
quantum black hole. 




(7.36) 
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